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SVANTE JANSON 

Abstract. Consider an election where A'^ seats are distributed among 
parties with proportions pi, ... , pm of the votes. We study, for the com- 
mon divisor and quota methods, the asymptotic distribution, and in 
particular the mean, of the seat excess of a party, i.e. the difference be- 
tween the number of seats given to the party and the (real) number Npi 
that yields exact proportionality. Our approach is to keep pi , . . . , pm 
fixed and let A — >■ cx), with A random in a suitable way. 

In particular, we give formulas showing the bias favouring large or 
small parties for the different election methods. 



1. Introduction 

The idea of a proportional election method is that each party gets a num- 
ber of seats that is proportional to the number of votes. The same math- 
ematical problem arises if seats are to be apportioned (before the election) 
between states in a union or between multi-member constituences according 
to their populations. In particular, the problem of apportionment to the 
House of Representatives in the United States has been the source of much 
debate as well as much research since 1790, including the invention of sev- 
eral important election methods; see Balinski and Young for a detailed 
history. (The European Union has yet to agree on a similar procedure; see 
Grimmett flo| for a recent attempt to replace the current political bartering 
with a "mathematical formula", i.e. a well-defined method to allocate seats 
in the European parliament to the member states.) To fix the terminology, 
we talk in this paper about parties in an election; the case of apportionment 
between states etc. differs mathematically only in language. 

Of course, exact proportionality is in general not possible, since the num- 
ber of seats has to be an integer, and the election method can be seen as 
a method to "round" the real numbers that yield exact proportionality to 
integers. The purpose of this paper is to study the resulting deviation from 
perfect proportionality, i.e. the "rounding errors" introduced by the method. 
In particular, we are interested in whether there is a systematic bias favour- 
ing large or small parties for various election methods. 

This problem has received a considerable amount of attention going back 



to (at least) Polya |23l. |24|. |25|. (Sainte-Lague [27j did a related comparison 
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between two different metliods already in 1910, see Example 13.51 ) Some in- 
vestigations have analysed data from real elections or from simulations (with 
party sizes chosen at random according to some predetermined distribu- 
tions), see e.g. Balinski and Young [1] and Schuster, Pukelsheim, Drton and 
Draper |28| . The problem has also been studied theoretically by assuming 
that the parties have random sizes with a uniform distribution of the relative 
sizes over all possibilities, see Polya (23, 24, 25|, Schuster, Pukelsheim, Drton 
and Draper j28|, Heinrich, Pukelsheim and Schwingenschlogl jl3i] and Drton 
and Schwingenschlogl (sO, 0, |3l|. Random party sizes with more general 
distributions have been considered by Heinrich, Pukelsheim and Schwingen- 
schlogl ||12|. [l3| and Schwingenschlogl {2^. (See also Section |6] below. ) 

The approach in the present paper is different. We consider a number of 
parties with fixed sizes and let the total number of seats N be random. More 
precisely, we let be chosen uniformly at random in 1,2,..., Nq for some 
large integer A^O) and then we take the limit as Nq^ oo. (Limits obtained 
in this way, if they exist, can be regarded as results for "a random positive 



integer" N. See also Remark 13.11 ) The same approach has been used in [17 1 
for a related problem (the Alabama paradox for Hamilton's method). 

We use this approach for several different election methods, and calcu- 
late, for given sizes of the parties (under a technical condition, see Sections 
|3] and [TT]) . the (asymptotic) distribution of the deviation from perfect pro- 
portionality as well as its mean, i.e. the systematic bias, and its variance. 
In particular, this shows how the bias depends on the sizes of the parties in 
a more detailed way than previous studies that consider random party sizes 
(see the references given above). 

The election methods that we consider are divisor methods of the "linear" 
(or "stationary") type (including Jefferson/D'Hondt and Webster /Sainte- 
Lague) and quota methods (including Hamilton/Hare and Droop). See fur- 
ther Section [21 and Appendices IAVIbI 

The main results are stated in Section |3] and proven in Sections HHSl Sec- 
tion |6] show that the results hold also in the more traditional approach with 
deterministic house size A^ and random party sizes. The following sections 
contain some simple applications of our main results: Section [7] discusses 
the probability of violating quota. Section [8] the expected gain or loss for 
parties forming an apparentement (coalition) , Section |9] the Sainte-Lague 
divergence studied by Heinrich, Pukelsheim and Schwingenschlogl jlH, Il3l | 
and Section [10] some related goodness-of-fit functionals. Section [TT] discusses 
our technical condition and the case of rational sizes of the parties. The 
appendices contain some background material. 



Acknowledgement. I thank Friedrich Pukelsheim for many helpful com- 
ments. 
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2. Notation 

We assume throughout that we have m parties with Vi votes for party 
i; we let V := YlT^i'^i the total number of votes and pi := Vi/V the 
proportions of votes for party i. (In real elections there may also be votes 
for parties that are disqualified because of a threshold rule, blank votes, and 
other invalid votes; these are ignored here so Yll^iPi ~ consider only 

parties that have received at least one vote, so pi > 0.) We further assume 
that seats are to be distributed (the house size), and let Sj be the number 
of seats given to party i; thus X^^^x Si = N. We occasionally write Si{N) 
when we need to specify the house size. 

Strict proportionality would give 

qi:=^N=piN (2.1) 

seats to party i. (This is usually not an integer.) We define the seat excess 
for party i to be the difference 

A,:=Si-qi = Si-piN. (2.2) 

Note that 

m m 

5^A, = ^Si-7V = 0. (2.3) 

i=l i=l 

We use the standard notations \x\ and [x] for rounding down and up of 
a real number x, i.e. the largest integer < x and the smallest integer > x, 
respectively. We denote the fractional part of x by {x} := x — [xj. 

More generally, let a be a real number. We say that the a-rounding of a 
real number x is the integer [x]^ such that 

X - a < [x]^ < X - a + 1; (2.4) 

equivalently, 

N„ + a-l <x< [x]„ + Q. (2.5) 

If X — a is an integer, we regard [x]^ as multiple-valued, and accept both 
X — a and x — a + 1 as values of [x]^. (This exceptional case typically occurs 
at ties in the election methods.) Consequently, 

[x]^= [x-a] = La: + l-aJ, (2.6) 

except the case when x — q is an integer and both x — a and x — a + 1 
are possible values. If < a < 1, this means that x is rounded down if its 
fractional part is less than a and up if its fractional part is greater than a. 
In particular, a = \ yields standard rounding, a = yields rounding up and 
a = 1 yields rounding down. (But note that we allow also a < or a > 1, 
in which case |x — [x]^ | may be greater than 1.) 

The two types of election methods that we consider can be defined as 
follows; see further Appendices IAVIbI where further details and equivalent 
definitions are given. Note that all methods that we consider are homoge- 
neous, i.e., the result depends only on the proportions pi of votes, but not 
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on the total number V of votes; hence we will mainly use pi rather than Vi 
in the sequel. 

• The 13-linear divisor method, or the divisor method with divisors 
d{n) = n + /? — 1 (where /3 is a real number): Let 



'Vi 




\Pi 1 




.D. 




[d'\ 





(2.7) 



where D (or D' = D /V) is chosen such that Xll^i ^« = N . (If /3 < 
or /? > 1, there are some exceptions for small A^, see Appendix lA.il ) 
This includes several important methods, in particular /3 = 1 (Jef- 
ferson, D'Hondt), /3 = 1/2 (Webster, Sainte-Lague) , /3 = (Adams). 
• The j-quota method (where 7 is a real number): Let Q := V/{N + j) 
and let 

= m + J)p^U, (2.8) 

a 

where a is chosen such that X^i^i ~ principle, 7 can be any 

real number (provided A^ + 7 > 0); in practice only the cases 7 = 
(Hamilton, Hare, method of largest remainder) and 7 = 1 (Droop) 
are important, although also 7 = 2 has some limited use. 

Remark 2.1. In both divisor methods and quota methods ties may occur; 
then the a-rounding in (12. 7p or (12. Sp has two possible values (because x — aG 
TL in (12. 4p ) for (at least) two parties, with (at least) one of them rounded to 
the larger value and (at least) one rounded to the smaller value. We assume 
that ties are resolved at random, but ties will usually not be a problem for 
us because of our assumptions below. 

We let U(a, 6) denote the uniform distribution on the interval (a, &). We 
will use Ui and IJi for random variables with Ui ~ U(0, 1) and Vi ~ U( — ^, ^), 
and we assume that such variables with different indices are independent. 
Recall that then E C/^ = 1/2, E C/^ = and Var Vi = Var Vi = 1/12. 

We denote convergence in distribution by — convergence in probability 
by — ^ and equality in distribution by =. 

3. Main results 

As explained in the introduction, we assume that pi, . . . ,Pm are fixed, 
and let the total number of seats (house size) A^ be random and uniformly 
distributed on {1, . . . , A^o} for some large integer Nq\ we then consider limits 
as A'^o — )• 00. For convenience, we denote this situation by A^ 00. (Note 
that this implies A^ — ^ 00.) 

Remark 3.1. We assume that A^ is distributed in this way for simplicity, 
but the results can easily be extended to more general sequence of random 
A^ G N. For example, we may take A^ uniformly distributed on an interval 
[A^i(A;), A^2(fe)] and let /c — )■ 00, for any two sequences Ni{k) and N2{k) of 
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positive integers with N2{k) — Ni{k) — )• oo. In fact, we conjecture that the 
results hold for any that satisfies the weak condition 



ip{e) ■.= Ee'^'"^ ^0 for any 61 G (0,27r). (3.1) 

(This implies N — t- oo and excludes parity restrictions.) Note that Weyl's 
Lemma r4. ll extends to this situation. However, we have not verified all details 
for this case, and we leave further investigations to the interested reader. 

In order to obtain simple asymptotic results, we will usually make one 
mathematical simplification: we assume that pi,...,pm are linearly inde- 
pendent over Q, i.e., that there is no relation 

aipi H h amPm = (3.2) 

with all coefficients Oj rational and not all 0. (Equivalently, there is no 
relation (|3.2p with integer coefficients Oj, not all 0.) The case when this 
assumption does not hold is discussed in Section [11] 

Remark 3.2. Mathematically, this assumption is reasonable, since if we 
choose pi,...,pm at random (uniformly given that their sum is 1, as in 
Section |6]), they will almost surely be linearly independent over Q. However, 
for real elections the assumption is clearly unreasonable since vote counts Vi 
are integers and the pi thus rational numbers. Nevertheless, the results below 
are good approximations if the numbers pi have large denominators and there 
are no relations (13. 2p with small integers ai, . . . ,am, see Section [TT] For 
practical purposes we thus can use the results below as good approximations 
for any pi and any N that is large enough. (However, we do not investigate 
the rate of convergence, and how large N has to be.) 

3.1. Divisor methods. We begin with a simple deterministic bound, valid 
for all house sizes and numbers of votes, cf. [ISj, Satz 6.2.11]. It follows from 
the proof, or from Theorem 13.71 below, that the bounds are best possible. 
Proofs of the results below are given in Section |4] 

Theorem 3.3. Consider the ^-linear divisor method. Then 

Pi - 1 + (/? - 1) {mpi - 1) < Ai < (/? - 1) {mpi - 1) + (m - l)pi. (3.3) 

Equivalently, 



/^i-[P-\]{mp,-l) 



1 m-2 , , 

< 2 + ^P^- (3-4) 



For the asymptotic bias, we have a simple formula. 



Theorem 3.4. Consider the /3-linear divisor method, and suppose that pi, 
. . . ,pm are linearly independent over Q. Then, as N oo, 

EAi^ {(3-^){mpi-l). (3.5) 
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Note that the asymptotic bias in (13. Sp for a party depends only on its size 
Pi and the number of parties, but not on the distribution of s izes between 



the other parties. The hmit in (13. 5p agrees with the formula by iKopfermann 
[lil . (S) on p. 208] for the advantage (or disadvantage) for different parties, 
derived by geometric considerations. 

In particular. Theorem 13.41 shows that the bias is for every party when 
(5 = 1/2 ( Webster /Sainte-Lague). It is well-known that this divisor method 
is (asymptotically) unbiased, see e.g. for various justifications 

of this; it is satisfying that our approach confirms this, and shows that the 
method really is unbiased for a party of any size. 

For j3 > 1/2 (for example Jefferson/ D'Hondt with (3 = 1), the bias is 
positive for pi > 1/ra and negative for pi < 1/m. It is well-known that 
these methods favour larger parties; we here see exactly how much parties of 
different sizes are favoured, and that the boundary is exactly at the average 
size 1/m. Note also that the bias grows with m, especially for large parties. 
For example, for /? = 1, a large party with pi ~ 1/2 has bias ~ (m — 2)/4, 
while a small party {pi <^ 1/m) always has a bias ~ —1/2. 

For 13 < 1/2 (for example Adams with /? = 0), we have the opposite biases; 
the boundary is still at pi = 1/m. 



Example 3.5. Sainte-Lague [27| compared in 1910 the results of his method 
(which he called the method of least squares) and D'Hondt's method in the 
case of two parties (that he thought of as two of several parties in an election) 
with pi > P2. He did not consider the biases in our sense, comparing the 
distribution of the seats to the distribution of the votes, but his approach 
means studying the difference of the two biases. (Since we know that Sainte- 
Lague's method (/3 = 1/2) is unbiased, this is equivalent to calculating the 
bias of D'Hondt's method (/? = 1).) 

Sainte-Lague found that if k := P2/P1 G (0, 1], then the largest party gets 
on the average (1 — A;)/(2(l -|- k)) seats more with D'Hondt's method. This 
agrees with Theorem 13. 4[ with m = 2 and pi = 1/(1 -\- k). 

Sainte-Lague ||27,] also considered the case of random pi, assuming that 
the ratio P2/P1 is uniformly distributed over (0, 1], and found by integration 
the average gain 

for the larger party with D'Hondt's method. (Note that this differs from 
the average bias 1/4 = 0.25 for the largest of two parties derived by |28| 
assuming that p2 is uniformly distributed on (0, 1/2), see Theorem l6.4l below 
with m = 2. The reason is that the two probability distributions of the party 
sizes are different; they are uniform in different ways.) 



Remark 3.6. Sainte-Lague f27|| seems to have thought that this calculation 
for two parties applies to any two of several parties, and that therefore, if 
party A is larger than B, then B would on the average gain 0.193 seats 
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compared to A if the method is changed from D'Hondt's to Sainte-Lague's, 
regardless of the number and sizes of the other parties. Theorem 13.41 shows 
that this is incorrect; in fact, if A and B both are small parties, they both 
have essentially the same bias —1/2 with D'Hondt's method. It is true that 
for divisor methods, the distribution of seats to two parties A and B is 
the same as if the total number of seats for these two parties is distributed 
between them by the same divisor method, see Theorem IA.19I However, if 
both A and B are small parties, they will by (13. Sp get on the average almost 
1 seat less together than their proportional share, and if their total number 
is redistributed between them, most of this loss will be borne by the larger 
party, which offsets the advantage of the larger party in the redistribution 
and explains the fallacy. More generally, parties 1 and 2 will by (|3.5p on the 
average get together A^' := N{pi +P2) + (/3 — 1/2) (mpi + ■mp2 — 2) seats. If 
we distribute this number of seats between the two parties, the same formula 
would yield that party 1 on the average gets, with p'^ := pi/{pi +P2), 

N'p'i + (/3 - i)(2p; - 1) = iVpi + (/3 - i) (p'i(m(pi + p^) - 2) + 2p'^ - l) 

= Npi + (/3 - \){mpi - 1) 

seats, in accordance with Theorem 13.41 

We can describe asymptotically not only the mean but also the distribu- 
tion of the seat excess. We also obtain a joint limit distribution, which is 
explicit although a little involved. 

Theorem 3.7. Consider the f3-linear divisor method, and suppose thatpi, ■ ■ ■ ,P; 
are linearly independent over Q. Then, as N 00, for each i, 

m-2 

A, A X, := i/3 - l){mpi -l) + Uo+piY, Uk, (3.7) 

fc=i 

where ~ U(— ^, ^) are independent. Moreover, jointly for i = 1, . . . ,m, 

A, A X„ (3.8) 

where the limit random variables Xi can be constructed as follows: Let 
Ui, . . . , Um ~ U(0, 1) and let J £ {1, . . . , m} be random with P( J = j) = pj, 
j = 1, . . . ,m, with all variables independent. Let 

Vr.= 1{J ^ ^}U, = l]!^'' (3.9) 
[0, i = J, 

and let, finally, 

m 

Xi := p^ ^3 -Vi + {P- l){mp, - 1). (3.10) 
3=1 

Of course, as we shall verify in Section HI Xi = Xi. 
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The description in (l3.9p -( f3?T0]) of the joint distribution of Xi, . . . , 
is exphcit but rather involved, and it is quite possible that there exists a 
simpler description; cf. Remark 110.31 For m = 2, X2 = —Xi, and thus 

{Xi,X2) = (Xi, -Xi) where Xi := (2/3 - l)pi -/3 + Uohy We leave 

it as an open problem to try to find alternative expressions for m > 3. 

Corollary 3.8. Consider the ^-linear divisor method, and suppose that 
pi, . . . ,pm are linearly independent over Q. Then, as N — > 00, 

VarA.^i±i!^i^ (3.11) 
12 ^ ^ 



and, when i 7^ j. 



Cov(Ai,Aj)^ (3.12) 

Note that the asymptotic distribution of the seat excess depends on f3 
only through its mean EXj (i.e., the asymptotic bias in (13. Sp ): the centred 
random variable Xj — EXj does not depend on /3. We see also, from (13.10p or 
(j3.1ip , that this random part of the seat excess tends to be larger in absolute 
value for a large party. 

In the limit as pi — )• 0, i.e. for very small parties, (|3.9p and (|3.10p yield 

Vi Ui and 

Xi^-Ui-i/3-l) = l-Ui-f3 = Ui-(3^ U(-/3, 1 - f3). 

Remark 3.9. A simple heuristic motivation for the bias in Theorem 13.41 
is that the /3-rounding in (|2.7p is unbiased if /3 = 1/2, but otherwise, it 
costs on the average each party /3 — 1/2 seats, compared to Vi/D which is 
truly proportional. In total, this makes a loss of m(/3 — 1/2) seats, which are 
redistributed proportionally (by shifting D), so party i gets back pim{f3— 1/2) 
seats, on the average, giving a net gain of {mpi — l)(/3 — 1/2), as shown by 
Theorem 13.41 

For a related, but rigorous, argument, let us compare the methods with 
parameters /3 and /3 + 1, temporarily using a superscript to distinguish the 
methods. (For example, Adams's method and Jefferson's method.) By The- 
orem IA.12I (assuming that is large enough) , 

sf\N) = sf^'^\N-m) + l. (3.13) 

Hence, 

Af ^ ( A^) = sf ) (N) - Npi = Af {N -m)-mpi + l. (3.14) 

Taking the limit as A^ 00, we see that the asymptotic distributions also 
differ by mpi — 1. Hence, it is not surprising that the bias in Theorem 13.41 is 
mpi — 1 times a linear function of /3. 
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3.2.^Quota methods. We again begin with a simple deterministic bound, 



cf. |l8l . Satz 6.2.3]. It follows from the proof, or from Theorem 13.131 below, 
that the bounds are best possible. Proofs of the results below are given in 
Section |5l 

Theorem 3.10. For the j-quota method, 

/ l\m — 1 . / l\m— 1 

lipi--) <\<l(p^--)+ . 3.15 

V m/ m \ m/ m 

Again we have a simple formula for the asymptotic bias. 

Theorem 3.11. Consider the ^-quota method and suppose that pi, ■ ■ ■ ,Pm 
are linearly independent over Q. Then, as N oo, 

EAi^j(pi--). (3.16) 

As for the divisor methods, the asymptotic bias in (I3.16P for a party 
depends only on its size pi and the number of parties, but not on the dis- 
tribution of size s between th e other parties. The limit in (I3.16P agrees with 
the formula by KopfermannI fisl . (S) on p. 199] for the (dis) advantage for 



different parties, derived by geometric considerations. 

In particular. Theorem 13.111 shows that the bias is for every party when 
7 = (method of largest remainder/Hamilton/Hare). Again it is well-known 
that this quota method is (asymptotically) unbiased, see e.g. 



our approach confirms this and shows that the method is unbiased for a party 
of any size. 

For 7 = 1 (Droop), the bias is positive for pi > 1/m and negative for 
Pi < 1/m, just as for Jefferson's method discussed above; hence Droop's 



method too favours larger parties. (See |l8l . p. 118] for figures illustrating 
the cases m = 2 and m = 3.) Note, however, that the bias for Droop's 
method does not grow with m. A comparison between (|3.16p (with 7=1) 
and (|3.5p (with /3 = 1) shows that the bias for Jefferson's method is m/2 
times the bias for Droop's method, for any party of any size. Hence Droop's 
method is less biased than Jefferson's for every m > 3. (For m = 2 not only 
the biases are the same; in fact the methods coincide when there are only 
two parties, see Appendix IB. 11 ) 

Remark 3.12. There is a simple heuristic explanation of the formula (|3.16p 
for the bias, if we accept the fact that Hamilton/Hare's method (7 = 0) 
is unbiased. Consider 7 = 1 (Droop); this method can be seen as doing 
Hamilton's method with -|- 1 seats, but retracting the last seat. Since 
Hamilton's method is unbiased, this first distributes on the average [N -\-l)pi 
seats to party i. The retracted seat is taken essentially uniformly at random, 
since it depends on the fractional parts of the numbers qi = Npi only; hence 
each party loses on the average 1/m seats in that step. This argument gives 
a bias of — 1/m, in accordance with (I3.16p . 

Again we have an explicit expression also for the asymptotic distribution 
of the seat excess. 
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Theorem 3.13. Consider the j-quota method and suppose that pi, 
are linearly independent over Q. Then, as N oo, for each i, 



Ai^Yi:=j(p,--)+Uo + -\2Uk, (3.17) 



k=l 

where ~ U(— ^, ^) are independent. 

Corollary 3.14. Consider the 'j-quota method and suppose that pi, . . . ,pm 
are linearly independent over Q. Then, as N — > oo, 

l + (m-2)/m^ (m + 2)(m-l) 
Vai A, ^ = . (3.18) 

As for the divisor methods in Theorem 13.71 the asymptotic distribution 
of the seat excess depends on 7 only through the mean in (|3.16p ; 1^ — E 1^ 
does not depend on 7. Moreover, unhke the case of divisor methods, this 
centred random variable does not depend on the party size pi. In particular, 
for Hamilton/Hare's method (7 = 0), the asymptotic distribution of the seat 
excess is the same for every party. 

Remark 3.15. Let for simphcity 7 = 0. Note that the hmit variable Yi in 
(|3.17p is uniform C/(— ^, ^) for m = 2; for m > 3 the distribution is more 
complicated, but as m — >• 00, the limit converges to U{—^,^) (by (|3.17p 
and the law of large numbers). Hence, for large m, Yi has essentially the 
same distribution as for m = 2, and although the supremum and infimum 
of the range of Yi tend to 1 and —1 as m — t- 00, cf. (|3.15p . the probability 
P(l-^[-i,i])^0. 

Similarly, it is seen from the variance formula (|3.18p that the variance is 
1/12 PS 0.08333 for m = 2, increases to 5/54 0.09259 for m = 3 and to the 
maximum 3/32 = 0.09375 for m = 4; the variance then decreases again and 
converges to 1/12 as m — )• 00. 

The difference from a uniform variable is not very large, and it may in 
practice be a good approximation to regard the seat excess Aj in Hamilton's 
method as U(— ^, ^) for any party and any reasonably large house size. 

We have stated Theorem 13.131 onl v for individual parties. It is easy to see 
that a joint asymptotic distribution exists, but (unlike the case of divisor 
methods, cf. Theorem 13. 7p we do not know any simple description of it. 
(There is a description implicit in the proof of Theorem 13.161 see (|5.14p .) 

Theorem 3.16. Consider the 'j-quota method and suppose that pi, . . . ,pm 

are linearly independent over Q. Then, as N — > 00, Aj — > Yi jointly 
for i = 1, . . . ,m for some random variables Yi such that the distribution of 
(Yi — 7pi, . . . ,Ym — 7Pm) "is symmetric (exchangeable) and independent of 

pi, . . . ,Pm, md Yi = Yi for each i. 
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Corollary 3.17. Consider the 'j-quota method and suppose that pi, . . . ,pm 
are linearly independent over Q. Then, as N oo, when i ^ j, 

Cov(A„ A,) ^ Cov(y„ Y,) = --—1 = (3.19) 

Problem 3.18. Find an explicit formula for the limit variables Yi, . . . ,Ym 
in Theorem 13.161 



4. Proofs for divisor methods 



Proof of Theorems \3.3\ and \3.'/l For the /3-hnear divisor method, the seat 
distribution is given by (|2.7p with D and D' = D/V such that Yl^i — ^■ 
We let D' =l/t; thus, using (|2Z3) and HIM . 

Si = [Pit]fi=[Pit- P + (4.1) 

(with the usual possible exception if pit — /3 is an integer). Let here t grow 
from to oo; think of t as time. The seats are assigned one by one (unless 
there is a tie). 

Fix one party j and let tjr be the time when party j receives its r:th seat. 
By ()4T]) . r = pjtjr - (3 + 1, and thus 



t 



r + /3- 1 
Pj 



(4.2) 



At this time tjr, by (|4.ip again, party i has Sj seats with 



Si = [pitjr-(3 + l\ = 

Let, for i = 1, . . . ,m, 

Wi := {pitjr -13 + 1} 



^(r + /3-l)-/3 + l 
iPj 



Pi 



Pj 



P: 



(/3_1 



Pj 



Pi , (Pi 

—r + 



Pj 



Pj 



1 (/3-1) e[o,i] 



(4.3) 



(4.4) 



defining Wi := 1 in the exceptional case when Pitjr — ji G Z and Sj = pitj^ — ji. 
(But PVj := when Pitjr — /3 S Z and Si = Pitjr — /3 + 1.) Then (14. 3p can be 
written 



Si = ^r+(^-l){P-l)-Wi, i = l,...,m. 



Pj 



Pj 



(4.5) 



Note that Wj = 0, since Sj = r = Pjtjr — /? + 1. 

The total number of seats at the moment tjr is, by (14.5 



i=l 



Si = — r +1 m 

Pj ^Pj 



{f3-l)-J2Wi., 



(4.6) 



i=l 
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and thus the seat excess is, combining (14. Sp and (14. 6p . 

m 

^i = Si-Npi=(^-l){p-l)-W^-(^-mpMp-l)+piY,Wk 

\ Da / \ T)a / 



k=l 



Pj ' ^Pj 
m 

= P^Y,Wk-Wi + {|5-l){mp^-l). (4.7) 

k=l 

So far everything is deterministic. Theorem 13.31 follows from (14. 7p since each 
Wk £ [0, 1] and Wj = 0. 

Let us now assume that pi, . . . ,Pm are linearly independent over Q; in 
particular, this implies that Pi/pj is irrational when i ^ j- Note that there is 
a tie between parties i and j at time t if and only if both pit — (3 and pjt — f3 
are integers. If parties i and j tie at two different times ti and t2, we thus 
have Pi{ti — t2) G Z and — ^2) £ and taking the quotient we find 

Pi/Pj G Q) a contradiction. Hence there is at most one tie for each pair of 
parties, and thus at most (™) values of t (and thus of N) for which there is 
any tie. Since we consider asymptotics, we may ignore the ties and assume 
that no tie occurs. 

Let now be random in {1,... , iVo} and let J = J{N) be the party 
getting the A^:th seat. Then Aj is given by (14. 7p . where Wi = Wi{N) for 
each N is given by (|4.4p with j = J and r = sj{N). Fix again j and consider 
first only such that J = j. The number of such < Nq is Sj{NQ), the 
number of seats party j receives when the house size is A'^o; note that, e.g. 
by Theorem [331 Sj{No) = PjNq + 0(1) 



We recall a well-known result by Weyl |32||. (The result is usually stated 
with oi = • • • = Ofc = 0; this case implies immediately the more general 
version given here.) We say that an infinite sequence {vn)n>i £ [0, l)'^ is 
uniformly distributed if the empirical distributions n^^ SnLi converge to 
the uniform distribution as no —s- 00, where 6y^ denotes the Dirac measure. 
(Recall that this means that if ^4 C [0, 1)*^ with X{dA) = 0, then #{n < no : 
Vn S A}/nQ — )• X{A), where A is the usual Lebesgue measure, see e.g. jij.) 

Lemma 4.1 (Weyl). Suppose that yi, . . . , and 1 are linearly independent 
over Q, and let ai, . . . , be any real numbers. Then the sequence of vectors 
{{nyi + ai}, . . . , {ny^ + afc})„>i G [0, 1)'^ is uniformly distributed in [0, 1)'^. 

□ 



(The standard proof is by showing that the Fourier transform (charac- 
teristic function) ^ X]n=i ^^P(^^i Sj=i ^ii^'^i}) ~^ 0' ^0 00, for any 
fixed integers £1, . . . not all 0, see for example 0, Exercises 3.4.2-3].) 

Since pi, ■ ■ ■ ,Pm are linearly independent over Q, so are pi/pj, ■ ■ ■ ,Pm/Pj, 
i.e. {pi/pj : i ^ j}U {!}. Taking yi := pi/pj and a, := (y^ - l)(/3 - 1), we 
thus see from Lemma |4. II that (Wi)ij^j, given by (14. 4p for r = 1, . . . , sj^Nq), 
are uniformly distributed vectors in [0,1)™""^. In probabilistic notation, 
regarding Wi = Wi{N) as random variables, conditioning on J = j and 
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recalling Wj = when J = j, as Nq — )• oo and thus Sj{NQ)^ oo, 

{iwi,...,w^)\j= j) A {i{i + m)Zv (4-8) 

Since P(J = f) = Sj{NQ)/NQ — )• pj as Nq — )• c«, the random variable J has 
asymptotically the distribution given in Theorem 13. 7| and it follows from 
(|4.8p that, with Vi defined as in Theorem 13. 7| 

{Wi,...,Wm) ^{Vi,...,Vm). (4.9) 

Consequently, (|4.7p yields Aj — ^ Xi, jointly for all i = 1, . . . , m, with Xi 
defined by (l3T^ . 

To obtain the simpler form (|3.7p for an individual Aj, let 

m 

:= Xi-i^- Dimp, -l)=piY^ Vj - - imp, + i. (4.10) 

i=i 

Then, if J = i, 

m— 1 

-'^i = Pi ^ Uj - \mpi + 1 = ^ [7j + i - ipj, 

and if J 7^ i, 

m-2 

Xi=pi^^ Uj + (pi - l)Ui - \mpi + 1 = Pi ^ [7j + (pj - 1)Uq - \pi. 
Hence, if we define 

^ _\piUm-l + \-\pi, J = i, ^^^^^^ 
m-2 

Xi=PiY, Uj + z, 
J=l 

with Uj {1 < j < m — 2) and Z independent. By (|4.1ip . conditioning on 
J, we have the distributions (Z | J = i) ~ U(| — pj, ^) and (Z | J 7^ i) ~ 
U(— ^, ^ — Pi); furthermore P(J = i) = pi, and it follows that Z ~ U(— ^, ^), 

so Z = C/q and thus 

m-2 

X,=p,Y,Uj + Uo. 

i=i 

The definition (|4.10p now shows that Xi = Xi, with Xi defined in (|3.7p . □ 

Proof of Theorem \3.4\ and Corollary \3.8[ Since Aj is bounded, e.g. by The- 
orem 13.31 (or by the simpler Theorem IA.17P , the convergence in distribution 

Aj Xi in Theorem 13.71 implies convergence of all moments and mixed 



then 
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moments, see jTH, Theorem 5.5.9]; hence it suffices to calculate the expecta- 
tions, variances and covariances of Xi. 

We have EC/j = 1/2 and thus by (13. 9p . since J and Ui are independent, 

Ey, = P(J/i)EC/i = (l-p,)i (4.12) 
Consequently, ^YaLi = \ YJiLii^ - Pi) = \{'r^ - 1) and, by (I3.10p . 

m 

EX, = p,E^ V,- - E^i + (/3 - l){mp, - 1) 
i=i 

m — 1 1— Pi /„ 1\, 

= Pi — ^ ^ H (/? - l)(mpi - 1) = - -j{mpi - 1). 

This shows Theorem 13.41 

For the variance and covariances, we have for any i 

EV;'=¥iJ^i)EUf = il-pi)l 

and, when i 7^ j, by independence, 

E{ViVj) = E(1{J / i}l{J / j}U,Uj) = P(J ^ {i, j})EC/,E[/,- = ^-P^^-P^ , 
Consequently, using also (I4.12p . 



cov(y„y,: 



34 12 

1 - Pi - Pi 1 - Pi 1 - Pj _ PiPj _ ^PiPj 



4 2 2 4 12 

We obtain from this also 



Cov(y. E y,) = E Cov(y., y,) = - 5: ^ 

i=i i=i 
1 + 2pj 3pj 1-Pi 



12 12 12 ' 

mm m m ^ ^ 

j=i 2=1 i=i «=i 

Consequently, (|3.10p yields 

m m 

Var(Xi) = Var V,) - 2p, Cov (f^, F,) + Var(Fi) 

2^-1 1 - Pi 1 + 2pi - 3p2 (m - 2)pf + 1 

= P,- 2n j !- = 

12 ^ 12 12 12 

and, for i ^ j, 

m m m 

Coy{Xi, Xj) = pipj Var 14) - p, Cov(y,-, ^ y^) - p^ Co^(yu 

k=l k=l k=l 

+ CoY{Vi,Vj) 
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m - 1 pi{l - pj) + pj{l - Pi) SpiPj 

= Pipj 

12 12 12 

_ (m - 2)piPj -pi- Pj 
12 

Corollary 13.81 follows. □ 

5. Proofs for quota methods 

Proof of Theorems \3.1(A and \3.13\ The proof for quota methods uses argu- 
ments similar to the proof in |l7| , and the reader might wish to compare the 
versions. (Only the case 7 = is treated in jl7'|, but that is only a minor 
simplification. For completeness we repeat some results from |17].) For no- 
tational convenience, we show the results for party 1, i.e. we take i = 1 in 
the proofs. 

By definition, Sj = [{N + 7)pj]^ for an a that makes X]j=i ~ ^ • 
Explicitly, by (1231) . 



Sj + a-l<{N + -i)pj<Sj+a. (5.1) 
Thus, using party 1 as a benchmark, 

Sj- si-l<{N + -f){pj - Pi) < Sj - si + 1. (5.2) 
Since Aj — Ai = Sj — si — N{pj — pi), (|5.2p yields 

Aj- - Ai - 1 < j{pj - pi) < Aj - Ai + 1. (5.3) 
Thus, summing over all j ^ 1, recalling that YlJLi — ^^^'^ YlJLiPj — 1) 

— mAi — (m — 1) < 7(1 — mpi) < — mAi + (m — 1). (5-4) 

Hence, 

7P1-- <Ai<7Ui-- + , (5.5) 

V m/ m V m/ m 

which shows Theorem 13.101 

Suppose now that all pi — pj are irrational. A tie between parties i and j 
can occur only if {N + ^)pi — {N + j)pj £ Z. If this happens for two different 
house sizes A'^i and A''2 , then both ( A^i + 7) (pi — Pj) G Z and {N2 + 7) {pi — 
Pj) G Z, and thus by taking the difference (A'^i — N2){pi — Pj) € Z, which is 
impossible if pi — pj ^ Q. Hence, our assumption implies that there is a tie 
between parties i and j for at most one value of A^, and thus at most (™) 
values of for which there is any tie. Since we consider asymptotics, we 
can ignore these A^ and assume that there are no ties. Hence we can choose 
a such that there are strict inequalities in (|5.ip . and thus there are strict 
inequalities in (15. 2p . 

By (|5.2p (with strict inequalities) [{N +^){pj—pi)\ G {sj — si — 1, Sj — si}. 
Let 

Ij := Sj -si- [{N + ^){pj - pi)\ . (5.6) 
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Then Ij G {0, 1} and h = 0. We have 

Aj - Ai = Sj - si- N{pj - pi) 

= Sj - si-{N + j){pj - pi) + -f{pj - pi) 

= ij-m+ i){pj - pi)} + i{pj -pi). (5.7) 

Summing over j we obtain, again recalhng XljLi ~ 0) 

m m 

-mAi = ^Ij- Y^^iN + j){pj - pi)} + 7(1 - mpi). (5.8) 

j=i i=i 

Let L := X]J=i ^j- f|5.8p . we then have the formula 

1 1 / \ 

Ai = 7(pi--)+- 5;{(iV + 7)(P,-Pi)}-^ • (5.9) 



Note that L is an integer with < L < m — 1 and, by (15. 6p . 

mm m, 

j=l j=l j=l 

m 

= N -msi-Y,liN + 7){pj - pi)\ 

m 

= N -"^KN + j){pj - pi)\ (mod m). 

Let M.odm{x) := m{x/m}, the remainder when x is divided by m. We thus 
have 

L = Mod™ (n -f2l{N + 7)(P, - Pi)j) • (5.10) 

Let Uj := Pj—pi and := ^{pj—pi); thus (A^+7)(pj — Pi) = Nyj+aj. Let 
us now assume that pi, ■ ■ ■ ,Pm are hnearly independent over Q. It is easily 
seen that then y2, ■ ■ ■ , Um and pj = 1 also are linearly independent over 
Q. (This can be seen as a change of basis, using a non-singular integer 
matrix, in a vector space of dimension m over Q.) 

We need the following extension of Lemma 14.11 slightly generalizing 17 , 
Lemma 4.3]. The proof is given in Appendix ICl 

Lemma 5.1. Suppose that yi, . . . ,yk and 1 are linearly independent over Q, 
and let ai, . . . ,ak he any real numbers. 

Let ^„ = Modm(?^ — l^j=iL"'yj + '3^iJ) S {0, . . . , m — 1}. Then the sequence 
of vectors 

{{nyi + ai}, . . . , {nyt + Ofc}, 4) G [0, 1)^= x {0, . . . , m - 1} (5.11) 
is uniformly distributed in [0, l)'^' x {0, . . . , m — 1}. 
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Note that L is (I5.10p equals £jv in Lemma l5.1| with y.j and Oj as chosen 
above. Consequently, (15. 9p and Lemma [5.11 show that 

A,A,(p.-i)+iff:j/,-iy (5.12) 



where Uj ~ U(0, 1) and L is uniformly distributed on {0, . . . ,m — 1}, with 

all vs 
thus 



all variable independent. Moreover, U2 = l-U2 and L + [/a ~ U(0, m), and 



L-U2 = L + U2-l = mUo - 1 = mUo + (m - 2) /2. 



Hence, (|5.12p can be written 



1 1 / 

-) + - E^. 

Vi=3 



Ai A 7(pi - - ) + - I > - mC/o 



with C/j := Uj — ^ for j > 1, which by Uq = —Uq proves Theorem 13.131 □ 

Proof of Theorem \3.11\ and CoroUary \3.14\ Since Aj is bounded, e.g. by The- 

orem l3?T0| the convergence in distribution Aj — ^ Yi in Theorem 13.131 implies 
convergence of all moments and it suffices to calculate the expectation and 
variance of Yi. This is immediate since EC/^ = and YaiUk = 1/12; thus 

by dam, 

EYi=^(pi--), 
Var Yi = h 



12 12m2 ' 

Theorem 13.111 and Corollarv 13. 141 follow. □ 

Proof of Theorem \3.1(A Let Wi := {{N + 'y)pi}. For any j, and 
1 = Yl^iPi linearly independent over Q, and thus the vectors (Wj)j-^j 
are uniformly distributed in [0, 1)™-^! by Lemma [4.11 Moreover, 

m m 

^VFiEE^(iV + 7>i = iV + 7EE7 (modi), 

i=l i=l 

SO Wj =j - Y^i^j Wi (mod 1) and 

Consequently, {Wi, . . . , Wm) (Vi, . . . , Vm) where each Vi ~ U(0, 1), any 
m — 1 of them are independent, and any such set determines the last vari- 
able by Vj = {'J — 'Yli^j Obviously, the distribution of (Vi, . . . , Vm) is 
symmetric. 
We have 

Aj - 7Pi = Si-{N + -f)pi = Si- [{N + -f)pi\ - Wi. (5.13) 
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By the construction of the 7-quota method, the numbers Sj — [{N + 'j)pi\ 
depend only on {Wi, Wm), say Sj - [{N + 'y)pi\ = fi{Wi, Wm) for 
some function /j, and it follows from (I5.13P that, jointly for all i, 

A. - A := n{Vu ...,Vm)- F„ (5.14) 

where the distribution of {Zi,...,Zm) is symmetric. (The m functions /j 
differ only by permutations of the variables.) Now let Yi := Zi + jpi. Then 

Ai Yi, and thus Yi =Yihy comparison with Theorem 13. 131 □ 

Proof of Corollary \3.17[ By Theorem 13 . 1 6 1 and the fact that each Aj is boun- 
ded, Cov(Aj, Aj) — )• Cov(Yi,Yj). Moreover, by symmetry, Cov(l^,yj) is the 
same for all pairs with i 7^ j. We have X^^^Lx ^ = as a consequence 
of ([23]), and thus 

m 

= Cov(Yi, ^ y,) = Var(yi) + ^ Cov(l-, Y,); 

hence by symmetry Cov(l^,l^) = —Yar{Yi)/{m — 1) = —YsLr{Yi)/(m — 1) 
for i ^ j- The final equality follows from (|3.18p . □ 



6. Random party sizes 

As said in the introduction, there is a long tradition of studying the 
bias of election methods theoretically with a fixed house size and random 
party sizes. More precisely, in this approach one (usually) assumes that the 
vector {pi, . . . ,pm) is random and uniformly distributed over the simplex 
©m := {{Pi, ■ ■ ■ ,Pm) S (0, 1)"* : ^^iPi = 1} for some given m. One then 
orders the party sizes pi, . . . ,Pm as p^^ > • • • > P[m] and considers the seat 
excesses A[]^] , • • • , A[^] of the parties in the same order. (Thus, A[]^] is the seat 
excess of the largest party, and so on.) Results for the mean E A[^] of the seat 
excess for the A;:th largest party (either exact formulas for a fixed or asym- 



totic formulas as — t- 00) have been obtained by Polya |23l.l24l.l25|. Schuster, 



Pukelsheim, Drton and Draper [is'l, Heinrich, Pukelsheim and Schwingen 



schlogl |l3l | and Drton and Schwingenschlogl [30|, |6| ; results on the variance 
are given by Schwingenschlogl and Drton [31,]. (We follow these papers and 
order the parties in decreasing order, cf. Appendi x iDl ) 

Heinrich, Pukelsheim and Schwingenschlogl fl2l . Il3l | and Schwingenschlogl 
(29! have also, more generally, considered other distributions of pi, ■ ■ ■ ,Pm 
(with an absolutely continuous distribution on &m)- Note also that Sainte- 
Lague (2?! considered a different distribution for two parties, see Example l3.5l 
(with P2/P1 uniformly distributed, which is not the same as p2 uniform). 

We cannot obtain these results for a given A^ directly from our results 
for random A^, nor can we obtain our results for fixed pi, ■ ■ ■ ,Pm from these 
results. However, by combining the methods, we can translate one type of 
results to the other and thus see that our main results transfer to the case 
of random pi, ■ ■ ■ ,Pm- (Now with random pi in e.g. (13. 7p and (13.17p .) 
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Theorem 6.1. Let pi, . . . ,pm be random with an absolutely continuous dis- 
tribution on the simplex &m, oind let N ^ oo (with N deterministic). 

(i) For the f3-linear divisor method, the conclusions (l3.7p -( f3?T0]) of The- 
orem hold, with Ui and Ui independent of pi, . . . ,pm and with 

^{J = j \ Pl,---,Pm) = Pj- 

(ii) For the j-quota method, the conclusions of Theorems \3.13\ and \3.16\ 

hold; in particular (|3.17p holds, with ~ U(— ^, ^) independent of 
each other and of pi, . . . ,Pm- 

Proof. Consider first the /3-hnear divisor method. The proof of fl3l . Theorem 
3.1(i)] shows that as — t- oo, 

(Ai,...,A^) A (6.1) 
for some random variables Zj; with the notation in [l3| . 

Z, = -{U,+q-^+ sgn{D,)m,{D,)) + c{q - \)W,. (6.2) 

Here Wj is our pj, c is our m, q is our /?, Uj is our Uj and sgn{Dc)mj{Dc) is 
an explicit but rather complicated function oi Ui, . . . , Um and Wi, . . . , Wm] 
however, the point is that for us it suffices to know the existence of some vari- 
ables Zj such that (|6.ip holds; we can ignore the construction. The proof of 
(|6.ip is based on the fact that for any (deterministic) sequence f„ — )• oo, the 
vectors of fractional parts {{i^nPi}, ■ ■ ■ , {i^nPm-i}) are uniformly distributed 
in [0, l)'""^, see Lemma \C.1\ which is the analogue of Lemma 14.11 in the 
present context. 

Consider now uniformly random on {1, . . . , Nq} and let Nq — t- oo, i.e., 
suppose oo in the notation used in the previous sections. Then — ^ 

oo, so by conditioning on A^, it follows that (|6.ip holds for such random A^ 
as well. On the other hand, since the distribution oi pi, . . . ,pm is absolutely 
continuous, almost surely pi, . . . ,Pm are linearly independent over Q. Thus, 
by conditioning on pi , . . . , pm we may apply Theorem 13.71 which shows that 

(Ai,...,A^) A(Xi,...,X^), (6.3) 

where Xi, . . . , Xm are given by (|3.9p - (|3.10p (with our random pi, . . . ,Pm)- 
In the case of random pi, . . . ,Pm and random A^ oo, we thus have both 
(j6.ip and (|6.3p ; consequently, 

{Z^,...,Zm) = {Xi,...,Xm). (6.4) 

Consider again deterministic N ^ oo. We know that (16. ip holds, so (16. 4p 

implies that (16. 3p holds in this case too. Finally, (Xi \ Pi) = {Xi \ pi) by 

Theorem 13. 7| and thus Xi = ATj, so (13. 7p too holds in the present setting. 

For the 7-quota method, the proof of Theorem 13.161 applies to the present 
situation as well, provided we use Lemma IC.ll instead of Lemma 14.11 In 
particular, (15.14P holds, with Vi,..., Vm ~ U(0, 1) independent of pi, . . . ,pm- 
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Further. (Yi \ pi, . . . ,Pm) = | Pi, ■ ■ ■ ,Pm), and thus Yi = Yi in the case of 
random pi, . . . ,Pm too. The results follow. □ 

Remark 6.2. By considering pi, ■ ■ ■ ,Pm uniformly distributed on a small 
neighbourhood maxj \pi—p^ \ < e of some {p^, . . . ^ ^m, and then letting 
e — )• 0, we recover the limit distributions in Section |3l with a somewhat 
different interpretation as a double limit. (In this version pi, ■ ■ ■ ,Pm can be 
any numbers, including rationals.) 

In particular, we can apply Theorem 16.11 to the case of uniformly dis- 
tributed pi, . . . , pm taken in decreasing order as above, since then (p[x] , • • • , P[m] ) 
has a uniform distribution on the subset &m,> '■= {{Pi, ■ ■ ■ jPrn) G &m ■ Pi ^ 
■ ■ ■ ^ Pm} of &m- This gives immediately the following theorem. (We leave 
the case of joint distributions to the reader.) 

Theorem 6.3. Let pi, . . . ,Pm be random and uniformly distributed on &m, 
and let N — )■ cxd. 

(i) For the /^-linear divisor method, for each j < m, 

m-2 

Ay] A X, := (/? - \){mp^j] -l) + Uo+ P[j] ^k- (6.5) 

k=l 

(ii) For the 'j-quota method, for each j < m, 

^b] ^ ■■= 7(pb1 --)+Uo + -Y.Uk. (6.6) 

k=l 

Here Uj. ~ U(— ^, ^) are independent of each other and of pi, . . . ,Pm- D 

As before, this theorem implies moment convergence because Ay^ is bounded. 
In particular, for the bias in this approach we obtain the following re- 
sult. This was (except for the case 7 7^ in (ii)) conjectured by Schuster, 
Pukelsheim, Drton and Draper [28l| and proven by Drton and Schwingen- 
schlogl [6] and Heinrich, Pukelsheim and Schwingenschlogl [13]. 

Theorem 6.4. Let pi, . . . ,Pm. be random and uniformly distributed on &m, 
and let N ^ 00. 

(i) For the /3-linear divisor method, for each j < m, 

EA[,] ^ EX, = (/3 - i)(mEpy] - 1) = (/3 - i) / ^ - - 1 j . (6.7) 

(ii) For the 'j-quota method, for each j < m, 

E Ay] ^EY,= 7(Epy] - _) = 2 / ^ _ _ 1 j . (6.8) 
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Proof. The convergence to EXj or EYj- follows, as just said, from Theo- 
rem [6]3l The first equalities in (16. 7p and (16. Sp follow immediately from (16. Sp 
and (|6.6p since 'E.Uk = 0, cf. (|3.5p and (|3.16p . It remains only to compute 
Epy]. This is a standard formula, and for completeness we give a proof in 
Appendix ini see (iD^el) . □ 

We similarly obtain results for the variance. This was found for m = 2 
and 3 (except the case 7 7^ 0) by Schwingenschlogl and Drton [31]. 

Theorem 6.5. Let pi, . ■ . ,Pm be random and uniformly distributed on &rn, 
and let N — t- 00. 

(i) For the /3-linear divisor method, for each j < m, 
Var Ay] ^ Var Xj = 1 (l + (m - 2) Ep^ ,) + (/3 - \fm^ Var(p[,]) 



1 m — 2 1 m — 2 / 1 

~ 12 ^ 12m(m -M) ^ + i2m(m + 1) 1 ^ 7 

\ m + 1 ^-^ i'^ m + 1 \^-^ I J 

(ii) For the 'j-quota method, for each j < m, 
Var A[,] ^ Vary, = - 1 + — ^ +7' Var(p[,]) 



12 V 

(m + 2)(m-l) 2 
12m2 ^ ^ 



/ 1 1 1 /J^ 1 

\ m(m + 1) ^ w?(m -|- 1) \^ i 



(6.10) 



Proof. Again, by Theorem 16. 3| we only have to calculate Yar Xj and Var Yj-; 
the calculations are similar and we give the details only for Xj. We use the 
standard decomposition 

VarX,- = E(Var(X,- \ pi, . . . ,p^)) + Var(E(l, \ pi, . . . ,p^)). (6.11) 



The conditional variance and mean are obtained immediately from (|6.5p . 
which yields, cf. (|33|) and (fSTTT) . 

E{Xj \ p,,...,pm) = {P- l){mpy] - 1), (6.12) 
1 + (m - 2)p2 

Var(X, bi,...,p„) = ^. (6.13) 

The first equality in (|6.9p follows from (|6.1ip - (|6.13p . It remains only to 
find Ep^^^.j and Var(p[j]); these are well-known and we give a computation in 
Appendix ID] for completeness, see (ID.6p - (ID.8p . □ 
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Example 6.6. For m = 3 and Jefferson/D'Hondt's method we obtain from 
([67} with /3 = 1, E(A[i], A[2], A[3]) found akeady by 

Polya l2^. 

For m = 3 and Droop's method we obtain from (|6.8p with 7 = 1, 

IE(A[i], A[2], A[3]) ^ 

Example 6.7. For m = 3, we have for the /3-hnear divisor method 
(Var A[i] , Var Apj , Var Ajg] ) ^ 

V2592 ^72^'^ 2^ '2592 + 72^^ 2^ '2592 ^72^^ 2^ ^ 
and for the 7-quota method 

(VarA[i],VarA[2],VarA[3]) ^ 

/ 5 13 2 5 7 2 5 4 2\ 
V54 648^ ' 54 648^ ' 54 ^ 648^ J ' 

These were found by Schwingenschlogl and Drton [31] (where, however, the 
case 7 7^ is not treated). 

Asymptotic covariances can be computed in the same way, using (I3.12p . 
(f3l9]) and fiDlOl) . 

In the same way as in Theorems 16. 4f[63| we can obtain results for the mean 
and variance of Ay^ for random pi, ■ ■ ■ ,Pm with other distributions than the 
uniform one. In particular, this gives the formulas by Heinrich, Pukelsheim 
and Schwingenschlogl jl^ and Schwingenschlogl jl^ for the asymptotic bias 
in terms of Epy] = ^{pj \ Pi > ■ ■ ■ > Pm) (assuming that {pi, . . . ,Pm) has a 
symmetric distribution). 

For example, Heinrich, Pukelsheim and Schwingenschlogl (l3| consider, as 
a special case, the case of a threshold t (with < t < 1/m) and uniform 
distribution on pi > • • • > pm > t. In this case, one has (as in (l3| ) {pi — 

l/m,...,pm-l/m) = {l-mt){pl-l/m, ... ,p^-l/m), where {pt,-.-,Pm) 
has a uniform distribution on &m,>, and the mean and variance are easily 
computed. 

7. The probability of violating quota 

We say that a seat assignment Si satisfies lower quota if Sj > [qi\ and 
satisfies upper quota if Sj < \qi\ ; it satisfies quota if both hold. In terms of 
the seat excess Ai = Si — qi, the assignment satisfies lower [upper] quota if 
and only if Aj > —1 [Aj < 1], and it satisfies quota if and only if |Aj| < 1. 

As is well-known, Hamilton/Hare's method always satisfies quota, while 
Jefferson's and Droop's methods satisfy lower quota and Adams method 
satisfies upper quota, see e.g. jH or Theorems 13.31 and 13.101 above. It is also 
well-known that Webster /Sainte-Lague does not always satisfy quota, but 
that violations are unusual in practice. 
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Theorems 13. 71 and 13.131 enable us to calculate the (asymptotic) probabili- 
ties that quota is violated for various methods. We give two examples. 

Example 7.1. Consider Jefferson/D'Hondt's method. Bv Theorem 13. 3 1 with 
f3 = 1, Pi — 1 < Ai < {m — l)pi; hence lower quota is always satisfied 
(as is well-known jl|), but upper quota can be violated for a party with 
Pi > l/(m— 1). Consider the case m = 3, and suppose that pi > 1/2. (Then 
P2iP3 < 1/2, so quota is always satisfied for the other two parties; thus at 
most one party can violate quota.) By Theorem 13.71 with /3 = 1, letting 
U- — - — JJ- 

Ai A Xi = |(3pi - 1) + Uo+PiUi = 2pi - Uo-piUi; (7.1) 
hence the asymptotic probability of violating (upper) quota is 

(2m - 1)"^ 

F{Xi > 1) = F{Uo + PiUi < 2pi - 1) = ^ (7.2) 

2pi 

(The set of allowed (Uq, Ui) is a right-angled triangle with sides 2pi — 1 and 
(2pi - l)/pi.) 

If we let {pi,P2,P3) be random and uniformly distributed on 63, as in Sec- 
tional then pi has the density 2(1 — pi) and thus the asymptotic probability 
that party 1 violates (upper) quota is 

/■I (2v-l)'^ 

/ — ^•2(l-p)dp = ln2-2/3w 0.026. (7.3) 

Ji/2 2p 

Consequently, the asymptotic probability that some party violates (upper) 
quota is 

31n2 - 2 f« 0.079, (7.4) 



as found by Niemeyer and Wolf [21 1 



Example 7.2. Consider Jefferson/D'Hondt's method, and a party i with 
three times the average size: pi = 3/m. Then the bias EXj = 1, by Theo- 
rem|33]with /3 = 1. It follows by ([3?Z| and symmetry that F{Xi > 1) = 1/2, 
so the (asymptotic) probability that the party violates quota is 1/2. For a 
larger party, the probability is even greater. 

Example 7.3. The Swedish parliament contains after the general election 
in 2010 8 parties: two large with 30% of the votes each and 6 small with 
5-8% percent each. The seats are in principle distributed by Sainte-Lague's 
method. (We ignore here complications due to the division into 29 con- 
stituencies and the system with adjustment seats, which in 2010 did not give 
complete adjustment because the number of adjustment seats was insuffi- 
cient.) 

We note from Theorem 13.31 (with /3 = ^) that the small parties always sat- 
isfy quota. In fact, (|3.4p shows that for Webster /Sainte-Lague, only parties 
with pi>l/{m — 2) can violate quota. 
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For the large parties we have (approximatively) pi = 0.3, and thus by 
(HZJ (still with (3 = \) Xi = Uq + ^2Y,l=iUk- An integration (using 
Maple) yields P(Xi > 1) = P(Xi <-!)?« 0.00045. Hence, for each of the 
two large parties the (asymptotic) probability of violating quota is 0.0009. 

8. Apparentements 

Suppose that two or more parties decide to form a coalition in the elec- 
tion, so that their votes are counted together. (In some election systems, 
parties can register such a coalition, called an apparentement, and continue 
to have separate names and lists. Otherwise, the parties can always appear 
on election day as one party with a common umbrella name. In any case we 
assume that this does not attract or repel any voters, so the coalition gets 
exactly as many votes as the parties would have had separately, and that all 
other parties remain unaffected with the same number of votes as before.) 

It is well-known that for a /3-linear divisor method with /3 > 1 (e.g. Jeffer- 
son/D'Hondt's method) parties can never lose by forming an apparentement; 
they will get at least as many seats together as if they appear separately. 
When /3 < (e.g. Adams's method), the opposite is true, and parties will 
never gain by forming an apparentement (conversely, they may gain by splin- 
tering), and for < /3 < 1 (e.g. Webster /Sainte-Lague's method), they may 
both gain and lose. For D'Hondt's method, this effect is politically impor- 
tant in reality, especially in small constituencies, and apparentements are or 
have been a regular feature in many countries, see e.g. 

Leutgab and Pukelsheim [l9| give a detailed study of the resulting gains 
in a set of real elections, including examples, statistics and theoretical values 
assuming that the party sizes are random (as in Section |6]). We can now give 
similar theoretical results for given party sizes. For simplicity, we consider 
only the case of a single apparentement of two parties; larger apparentements 
can be treated similarly as well as several apparentements (their effects are, 
asymptotically at least, additive). 

Let the parties be i and j, and consider the expectation of the gain Sjj — 
(sj -|- Sj), where sij is the number of seats the parties get together as an 
apparentement . 

Theorem 8.1. Suppose that pi, . . . ,pm are linearly independent over Q. // 
two parties i and j form an apparentement in the election, the mean of their 
seat gain Sij — Si — Sj satisfies, as N — > oo, 

Eisij -Si- Sj) ^ (/3 - ^){l-pi-pj) (8.1) 
for the (i-linear divisor method, and 

E{sij-si-sj)^j( T)=^—t TT 

for the ^-quota method. 
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Proof. If Ajj is the seat excess for the apparentement, then 

Sij -Si- Sj = N{pi +pj) + Aij - {Npi + Ai) - {Npj + Aj) = Aij - Aj - A^-, 

so E(sjj — Si — Sj) = E Aij — E Aj — E Aj , which by Theorems 13.41 and 13.111 
converges to, recalling that Aij is the seat excess in a contest with m — 1 
parties, 

(/3 - i) (^{m - l){pi + pj) - 1 - {mpi - 1) - {mpj - 1)^ 

and 

for the /3-linear divisor method and 7-quota method, respectively. □ 

One sees in the same way that these expected gains (which are negative if 
/3 < 1/2 or 7 < 0) are balanced by a loss for each other party A; of (/3 — 
for the /3-linear divisor method and 7/(m(m — 1)) for the 7-quota method. 

We have seen that apparentements are favoured when /?>l/2or7>0, 
but the effect is rather small, especially for quota methods. (Nevertheless, as 
said above, the effect is large enough to be politically important in reality for 
D'Hondt's method.) Note that for the /3-linear divisor method, the gain in 
(18. ip does not depend on the number (or sizes) ot the other parties, while for 
the 7-quota method, the gain in (|8.2p is largest for m = 3 or 4 and decreases 
as 0{l/m) for large m. 

Example 8.2. For Jefferson/D'Hondt's method, two small parties that form 
an apparentement gain at most 0.5 seats together, and two large parties gain 
less. 

Example 8.3. For Droop's method, any two parties that form an appar- 
entement gain at most g ~ 0.167 seats together when there are 3 or 4 parties, 
0.15 seats when there are 5 parties, and less if there are more. 

Problem 8.4. What are, for the different methods, the asymptotic distri- 
butions of the seat gain Sij — Si — Sj for an apparentement? 

We have so far studied the gain of the two parties combined, but the 
parties are probably more interested in how the gain is split between them. 
Typically, the seats given to an apparentement are distributed, in a sub- 
apportionment., between the participating parties by the same election method 
as for the main distribution (the super- apportionment). It seems reasonable 
that the gain then, on the average, is split between the parties proportionally 
to their sizes, cf. Theorems lA.llI and IB. 31 However, we shall see that this 
holds only for the divisor methods. 

Consider first the /3-linear divisor method. Theorem 13.41 shows that, 
asymptotically as 00, the number of seats the apparentement gets 

is on the average 

N{pi +Pj) + (/3 - i)((m - l){p, +pj) - 1). (8.3) 
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If we use Theorem 13.41 again for the distribution of these seats between i and 
j, we obtain that i gets on the average 

(a^(p. + (/3 - i) ((m - - 1)) + (/3 - i) (2-^ 

^ Pi I Pj ^ Pi I Pj ' 

= Npi + {13- \){{m - l)pi + 1) (8.4) 

seats, which, comparing to Theorem 13.41 for an election without the appar- 
entement, means an average gain for party i of 

<''-5'(sfi--''0 = sf^(''-^)(i-f--fi). (S'^) 

which indeed is the proportional share Pi/{pi +Pj) of the joint gain in (|8.ip . 

There is, however, a gap in this argument. In our model, the seat assign- 
ment is a deterministic function of N both in the super-apportionment and 
in the sub-apportionment, and the calculation above using expectations for 
random N assumes that there is no hidden correlation, where the numbers 
of seats given to the apparentement are seat numbers that tend to favour 
one of the two parties in the sub- apportionment. It seems very unlikely that 
there will be such a correlation under our assumption that pi , . . . , p^ are 
linearly independent over Q, but we have so far no rigorous proof, so (18. Sp 
should only be regarded as a conjecture for the asymptotic average gain for 
party i. 

Problem 8.5. Verify rigorously that the expected gain for party i converges 
to the value in (18. 5 p as N 00. 

Nevertheless, if we tentatively use (18. Sp . and regard it as a valid approx- 
imation for real elections with finite N (and arbitrary pi, ■ ■ ■ ,Pm), we can 
draw some practical conclusions for D'Hondt's method (/3 = 1). If the par- 
ties i and j both are small, then (18. 5p shows that the gain for party i is 
almost Pi/{pi +Pj)- Thus most of the gain of the apparentement goes to the 
larger party, and a party has very little to gain by an apparentement with 
a party that is substantially larger. In fact, if party i can choose between 
several parties to form an apparentement, we see, perhaps surprisingly, that 
it gets the largest gain by choosing the smallest partner. (This is based on 
asymptotics as — )• 00. For a real election with finite N there is a trade-off 
since a really small partner will not help significantly; the advantage can 
hardly be larger than Npj, so presumably only parties with qj = Npj being 
at least 1/2 or so will be useful partners. It would be interesting to study 
such effects for finite A^, but that is beyond the scope of the present paper.) 

Consider now instead the 7-quota method. One important complication 
is that quota methods are not uniform, in the sense of Appendix I A. 31 Thus, 
even if the apparentement gets the same number of seats as the two parties 
would have got together if they had appeared separately in the election, the 
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sub-apportionment may give a different distribution of tliese seats between 
the parties. 

We calculate as above, with the same reservation that a rigorous verifica- 
tion still is lacking. Theorem 13. Ill shows that, asymptotically as N oo, 
the number of seats the apparentement gets is on the average 

N{p^+PJ)+l[{p^+Pj)-^^)■ (8.6) 

If we use Theorem 13.111 again for the distribution of these seats between i 
and j, we obtain that i gets on the average 



(^N{pi + Pj) + 7 ({pi + Pj) 



\\ Pi , f Pi 1 

+ 7 



m-lJJpi+pj ^Pi+Pj 2 

/ m — 2 Pi 1 \ / s 

= Np, + ^[p, + 8.7 

V m — 1 Pi + Pj IJ 

seats, which, compared to Theorem 13 . 1 1 1 for an election without the appar- 
entement, means an average gain for party i of 

/m — 2 Pi m — 2\ 
7 7-^ ^ • (8.8) 



— 1 Pi + Pj 2m 

This is quite different from the proportional share Pi/{pi + Pj) of the joint 
gain in (18. ip . In fact, we see that for 7 > (e.g. Droop's method), a party 
gains only if Pi/ {pi+Pj) > (m — 1) / (2m). In other words, a party will lose if it 
forms an apparentement with a party that is only a little larger. Typically, 
the redistribution within the apparentement by the sub-apportionment is 
more important than the collective gain. 

9. The Sainte-Lague divergence 



Sainte-Lague [261, l27| based his proposal of his method on the fact that it 



is the least squares method minimising the sum 

m Ar o m , Ar ^2 ™- a 2 

-^E^^^)^^E^^E|: (-) 

i=l i=l i=l 

multiplying this by V we obtain the equivalent quantity 

™ A2 

^■■=^^ = T.Tr' (9.2) 

i=l 

which is called the Sainte-Lague divergence by Heinrich, Pukelsheim and 



Schwingenschlogl |l2l . |13|, where the asymptotic distribution is studied un- 
der the assumption of random party sizes (see Section |6|). We obtain cor- 
responding result in our setting as corollaries of the results in Section |3l 
(Unfortunately, we do not get as explicit result for quota methods as for 
divisor methods because Theorem 13. 161 is less explicit.) The most interesting 
case is the ^-linear divisor method (Sainte-Lague's method) [lH, since then 
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S and S are the minima over all allocations, but it is also interesting to see 
how much larger S is for other methods )l3| . 



Theorem 9.1. Suppose thatpi, . . . ,pm are linearly independent over and 
let iV ^ oo. 

i) For the f5-linear divisor method, 



S 



«±A^-(f(V. + ^-l))^ (9.3) 



.=1 « i.i 



where Vi , . . . , Vm are as in Theorem \3_ 
(ii) For the 'j-quota method, 

5A5^,^:=V^, (9.4) 

where Yi, . . . , are as in Theorem \3.1(A 

Proof. For the /3-linear divisor method, we use Theorem 13.71 and obtain by 
the continuous mapping theorem {l. Theorem 5.1] and (|9.2p immediately 

5 Y.T=i Xf/Pi- We write fICT]) as Xi = PiA-Ai, where Ai := Vi + f3-l 
and A := ZT=liV^ + /? " 1) = EI^i Thus 

1=1 i=l 1=1 i=l i=l 

which is M . 

The result for the 7-quota method follows directly from Theorem 13. 161 □ 

The expression in (19. 3 p for S'/j^m may be compared to the somewhat differ- 
ent expression in jl2i. Il3] for the limit for fixed and random pi, . . . ,Pm] by 
the argument in Section |6l the formulas have to be equivalent (for random 

Pi,...,Pm)- 

For the expectation we similarly obtain: 



Theorem 9.2. Suppose thatpi, . . . ,pm are linearly independent over Q, and 
let N -^00. 

i) For the /3-linear divisor method, 



ES 



(ii) For the 'j-quota method, 

12m^ Pi \i—i P^ J 
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Proof. This can be shown from Theorem 19. Ij but it is easier to use 

J^EAf _^ VarA, + (EA,)^ 

and Theorem 13.41 Corollary 13. 8| Theorem 13.111 Corollary 13.141 □ 

Note that YliLi ^/Pi ^ with equality if and only if all pi are equal. 
We know that S always is minimised by the i-linear diyisor method (Sainte- 
Lague's method). We see from Theorem 19.21 that also asymptotically, KS 
is strictly larger for any /3-linear diyisor method or 7-quota method, except 
when all parties haye exactly the same size (a triyial case where all methods 
yield the same result; note that we haye excluded this case by our assump- 
tions). 

Furthermore, Theorems 19.1H9.21 show that S is typically of the order m? 
if all Pi are of roughly the same order, but if some pi is yery small, S can 



12, 1 



be substantially larger. Heinrich, Pukelsheim and Schwingenschlogl 
consider S and the limit Sp^m for random, uniformly distributed, pi, ■ ■ ■ ,Pmi 
and show that, as m — )• 00, 

Sp,Jrr? -{^ + {(3- 1)2) log m A 5*, (9.7) 

where S* is a certain 1-stable random yariable (depending on (3). We see that 
the logm term and the large tail of S* (we haye ES"* = 00) depend on the 
possibility of some yery small pi] the asymptotic behayiour for deterministic 
pi, . . . ,pm with, say, {mmpi)~^ = 0{m) is quite different. To see this better, 
we first replace Vi (which are dependent) by Ui (which are i.i.d.) in (19. 3p : the 
following lemma shows that for large m, the difference is small. (Cf. similar 
approximations in (lil. [l3l|.) 

Lemma 9.3. Let Ui := Ui + (3 — 1, so Ui are i.i.d. U(/3 — l,/3), and let 

m jj2 /jn^ \ 2 



■ , P- V- 1 

1=1 1=1 



Then E|5^^rrt ~ S'i^m\ ~ 0{m), uniformly in pi, . . . ,pm. 

Proof We haye \Ui+(3-l\ < |/3| + 1 = 0(1) and similarly \Vi+P-l\ = 0(1). 
Further, by (|3J|), F{Vi / C/,) = P(J = i) = pi. Hence, 

^ \(V, + ,3-ir-(U, + fi-,f\ ^ OM ^ 0(1). (9.9) 

Pi Pi 

Similarly, + /? " 1) " ET=iiUi + P - I) = -Uj = 0(1) and thus 

/ m \ 2 / m x2 

(j2iVi + /3-l)] - K]([/, + /3-l) =0(m). (9.10) 
The result follows by summing (19. 9p and (19.10p . □ 
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We can now state a normal limit theorem corresponding to (19. 7|) . but for 
deterministic pi, ■ ■ ■ ,Pm- 

Theorem 9.4. Suppose that for each m, we are given pi, ■ ■ ■ ,Pm such that 
m^/^minj^j — ;> oo as m ^ oo. Let Am := Y.^iPT'^' '■= YlT=iPi'^ '^"■'^ 

m 

Cm '■= '^{Pi^ - "^)^ = Bm- 2mAm + . 
i=l 

Then, as m ^ oo, 

Sp,m - {Am/12 + {Am - m^){f3 - 1/2)2) ^ ^ 

Proof. Note first that, e.g. by Jensen's inequality, Am > rn?-, Bm > mn^ and 

Cm>^- 

Write, for convenience, 6 := /3 — 1/2. We have Ui = Ui + b and KU = b. 
Thus, (|9.8p can be expanded as 

^/3,m = 2^ \mh + }_^U^\ 

™ 7T2 _i_ oATT / \ 2 

= X; ( ^ ^ - 2m6^.) + b^^^ - b2^2 _ ^ ^ \ _ 

The last term has expectation Var Xll^i = 0{m), and can thus be ignored. 
Define _ _ 

Pi 

Then _ 

EU? 1 

EZi = = 

Pi I2pi 

and, by a straightforward calculation, 

180p2 3 \pi 

The random variables Zi are independent and each is bounded with \Zi\ = 
0{l/pi + m) = o{m?/'^) = o{Bm'^). The result now follows by the central 
limit theorem, using either Lindeberg's or Lyapounov's condition, see e.g. 
[nl. Section 7.2]. □ 

Remark 9.5. Sainte-Lague [27] also studied the least squares functional 
calculated per seat and not per voter: 

1=1 i=l 

and found that this is minimised by the divisor method with d{n) = ^Jn{n — 1), 
later suggested by Huntington, see Appendix |A] and Remark lA.lOl The 
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asymptotics of this quantity is, after division by {V/N)'^, the same as for S, 
because Si/N ^ pi. 

10. Some other goodness-of-fit functionals 

Many election methods, including the ones treated here, can be charac- 
terised as minimising various functionals, see e.g. jH, jTs'] and ^d\. Theorems 
13.71 and 13.161 yield by the continuous mapping theorem asymptotic distribu- 
tions for many such functionals for the methods treated here, under our 
standard assumptions that pi, . . . ,Pm are linearly independent over Q and 
N oo, which we assume below. We give some examples, leaving others 
to the reader. Note, however, that our results for quota methods are less 
complete, cf. Problem 13.181 

10.1. Hamilton divergences. Hamilton/Hare's method, the 0-quota me- 
thod, is the unique method minimising maxj |Aj|. It also the method min- 
imising maxj Aj. Moreover, it is, for any convex function (/? : M — [0, oo] with 
(/^(O) = and ip{x) > for x 7^ 0, the unique method minimising Y2i 
as was shown by Polya [23|. 

For the /5-linear divisor method. Theorem 13.71 implies 

max |Aj| — ^ max \Xi\ (10-1) 

l<i<m l<i<m 

with Xi given by (|3.10p . However, we do not know any really simple descrip- 
tion of this limit variable. 

For the 7-quota method, including Hamilton/Hare's method and thus the 
minimum of maxj |Aj| over all allocations. Theorem 13.161 vields 

max |Aj| — ^ max |yi|, (10-2) 

l<i<m l<i<m 

but we have no explicit description of this limit variable. 
Similar results hold for maxj Aj. 

We obtain more complete results if we instead consider the least squares 
function al T^. ^i A?, which also yields Hamilton/Hare's method by Polya's 
theorem |24| . (This case was shown earlier by Sainte-Lague |27| . who attrib- 
uted the result to Zivy.) 

Theorem 10.1. With assumptions and notations as in Sectionl^i for the 
P-linear divisor method, 

m m 

Y,^j^Sl^:=Y^Xl (10.3) 

i=l 1=1 

m ™ 1 _i_ / Q^i 2 m 

E A,^ ^ E Sl^ = Y: ^ + (/3 - \f Y^^rnp, " (^O'^) 

i=l 4=1 i=l 
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i=l 



:io.5) 



and for the j-quota method, 



m m 



^A? A5«^:=^y,^ (10.6) 



i=l i=l 
m 



E A? ^ ESH,,. . + + - l)^ (10.7) 

i=l i=l 

Proof. By Theorems ITTI and [3l6l calculating EXf and EF^^ rpj^g^, 
rem 13. 4| Corollary 13. 8| Theorem 13.111 Corollary 13.141 as in the proof of 
Theorem [921 □ 

10.2. Jefferson and Adams divergences. Heinrich and Schwingenschlogl 



14l | consider the functionals 



:= max (^-n) = max — (10.8) 

l<i<m\pi J l<i<m 



and 



:= min (^-n\= min — . (10.9) 

l<i<m\pi J l<i<m p^ 

(We could as well use the maximum and minimum of si/vi — N/V .) Since 
YT=i^i = 0> we have 5^ > but < 0, and we will therefore use 

It is easy to see that S'^ and \S^\ are minimised, over all allocations of N 
seats, by the Jefferson and Adams methods, respectively fg^l, [l|. 

Theorem 13.71 yields the asymptotic distribution of S'^ and | \ for any j3- 
linear divisor method; we consider for simplicity only the optimal methods. 

Theorem 10.2. With assumptions and notations as in Section 0, for the 
optimal allocations, 

m—l 

S' A 5i := ^ [/„ (10.10) 
i=i 



and 



We have 



\S^\^\St\:=m-^V,- mm (10.11) 

^ — ' l<t<m Pi 



lKSi = E\Si\ = ^^. (10.12) 
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Proof. S"^ is minimised by Jefferson's method, so Theorem 13.71 with (3 = 1 
yields, noting that minVi/pi = by (13. 9p . 

b := max > max — = > Vj— mm — 

l<i<m Pi l<i<m Pi ■^-^ l<i<m Pi 

m jTi— 1 

i=i jVJ i=i 

Similarly, 15"^! is minimised by Adams's method, so Theorem 13.71 with 
/3 = yields 

icAi . Aj d . Xi Y^TA , . 1-Vi 
\b I := — mm > — mm — = — > Vj+m— mm . 

l<i<m Pi l<i<m pi ^—^ l<i<m pi 

We have E 5";^, = (m - l)/2 directly from (110. lUI) . 

For \S:^\, we first note that if B := mini<j<m 1/Pi, then for < x < B, 



mm 



l-Vi 



> x] = F(Vi < 1 — xpi, i = 1, . . . ,m) 
l<i<m Pi y 

m 

= '^Pj f{Vi < 1 - xpi, i = 1,. . . ,m \ J = j) 

m 

= Y,P3^{U, <l-xpi, i^j) 

m 

= ^PjYlc^ - xpi), 



while P^mini<j<m(l — Vi)/pi >xj=OfoTx>B. Hence, 
in mm = / rl mm > x ] ax 

l<i<m Pi Jq \l<i<m pi 



,.B rn 



i=l 



XPi) 



1. (10.13) 



x=0 



Furthermore, Ej^J^^Vj = J2fJi Uj = (m - l)/2, and (ITOT^ follows 
from (fTOlTT) . □ 

Remark 10.3. Note that 5^ in (jlO.lOp does not depend on pi, . . . ,Pm, but 
for 1 5*^ I, the situation is more complicated; the representation in (jlO.lip de- 
pends on pi, . . . ,Pm, but it is still possible that has the same distribution 
for all pi, . . . ,pm and that we just have found an unnecessarily complicated 
description of it. Note that ElS*^] does not depend on pi, . . . ,Pm- In the 

case m = 2, it is easy to see that \S^\ = U = S'l^ for every {pi,P2), but we 
do not know whether something similar is true for m > 3. 
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Problem 10.4. Is jS"^] = for every pi, . . . ,Pm? 



The asymptotic results by |14| . taking the limit of the limits as m — t- oo, 
follow easily in our setting too. 



Theorem 10.5. For any pi, . . . ,Pm depending on m, as m ^ oo, 

^/^(5i-m/2) AiV(0, j^), 

V^(|5^|-m/2) AiV(0,^). 

Proof. The result for is immediate from (jlO.lOp and the central limit 
theorem. 

For 1 5^ I we use (jlO.lip . We have Yl^=i{^j — Vj) = Uj and consequently, 
by (fTOTTT) and (ITOT^ . 



E 



K\-Y.(^-Uj) 



< E 



^ -V- 



l<i<m V,: 

i=i - - 

< E C/j + E min ^~ ' < 2, 

l<i<m 



and the result follows by the central limit theorem applied to ^j^x(l — 
U,). □ 

11. Rational pi 

In our main results we have assumed that pi, . . . , pm are linearly indepen- 
dent over Q, since this is necessary for Lemma [4. II We consider briefly what 
happens when this is not satisfied, i.e., when pi, ■ ■ ■ ,Pm are linearly depen- 
dent over Q. In particular, we are for obvious practical reasons interested in 
the case when pi, ■ ■ ■ ,Pm are rational. (See also the corresponding discussion 
in 0.) 

Note first that it may now happen that there are ties. We assume that any 
ties are resolved randomly; thus the quantities Si{N) and Ai{N) in general 
may be random variables, also for a fixed N. This is no real problem, how- 
ever, and when discussing the means we may replace them by their average 
over all solutions in case of a tie. 

Lemma 14.11 does not apply when pi, . . . , pm are linearly dependent over 
Q, but the proof of it sketched in Section |4] shows that the sequence {{nyi + 
ai}, . . . , {nyk + afc})n>i is uniformly distributed on a coset of a subgroup 
of [0, l)'^; more precisely, if we for simplicity take all = 0, the empirical 
distributions converge to the uniform probability measure //y on a subgroup, 
with Fourier coefficients given by 

I (J otherwise. 

The proofs of Theorems 13.71 and 13.161 now show, with minor modifications, 
that the seat excesses A, still converge jointly in distribution as N — > oo. 
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but the limit distribution is different, and depends on pi,...,pm tlirougli 
some measures of the type Oll.ip . In particular, the bias limEAj is well- 
defined, but we do not have an explicit formula for it. 

However, if we consider a sequence of party size distributions {pik, ■ ■ ■ , Pmk), 
k = 1,2, . . . , (with a fixed m), such that pik — )• pi for each i as k ^ oo and 
further, for every integer vector (ai, . . . ,am) 7^ 0, we have Yl^i o-iVik 7^ 
for all large k, then it follows from (lll.ip and an inspection of the proofs of 
Theorems 13.71 and 13.161 that the measures ^Uy that now appear in the proofs 
will converge to the uniform measure /i, and thus the limit variables Xi and 
Yi will converge to the corresponding variables given in Section [3] for the case 
when pi, . . . ,Pm are linearly independent over Q. 

This shows that the results above are good approximations also for linearly 
dependent pi, ■ ■ ■ ,Pm, as long as there is no linear relation (|3.2p with small 
integers a,. In particular, for rational pi, ■ ■ ■ ,Pm we typically have a good 
approximation unless some pi has a small denominator. 

In the remainder of this section, we suppose that pi, . . . ,Pm are rational, 
with a common denominator L. Theorems lA.llI and IB. 31 then show that 
the sequence of seat excesses Aj := Si{N) — Npi, N > 1, has period L. (If 
/? < 0, /? > 1 or 7 < we may have to except some small N.) This gives 
another (simpler) proof that the limit distribution of Aj as N oo exists 
in this case, and shows that it is the same as the distribution if we take 
uniformly in a period {K + 1, . . . , K + L} (for any K that is large enough). 
In particular, the asymptotic distribution is discrete, and therefore not the 
same as in Section |3l The asymptotic bias lim E A j is easily calculated by 
taking the average of Aj over a period. 

Example 11.1. Let m = 2 and take {pi,P2) = (|j 5)- Then the ^-linear di- 
visor method equals the 7-quota method with 7 = 2/3—1, see Appendix lB.il 
Assume, for simplicity, < /3 < 1, i.e., —1 < 7 < 1. A simple calculation 
shows that for = 1,2,3, we have Ai = -^,—1,0. Hence the average is 0, 
so the methods are unbiased in this case, for any /3 E (0, 1) or 7 E (—1, 1)- 
(In fact, for this example the methods all coincide with Hamilton/Hare — 
Webster/Sainte-Lague.) In particular. Theorems 13.41 and 13.111 do not hold 
for rational pi, . . . ,Pm and (3 E (0, i) U (i, 1) and 7 E (-1,0) U (0, 1). 

If we instead take (3 = 1, i.e. 7 = 1, we have a tie for N = 2, with Ai = — | 
or |. The average of Ai over a period now is g. 

For the methods with (3 = 1/2 and 7 = 0, the result in Example 111.11 
agrees with Theorems 13.41 and 13.111 In fact, this is true for any rational 

Pl,...,Pm- 

Theorem 11.2. The ^-linear divisor method (Webster/Sainte-Lague) and 
the 0-quota method (largest remainder/Hamilton/Hare) are asymptotically 
unbiased for any rational pi, ■ ■ ■ ,Pm- 

Proof. Let, as above, L be a common denominator of pi, . . . ,Pm- It is for 
both methods easily seen, arguing similarly to the proofs of Theorems lA.llI 
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andlEl that Si{L-N) = Lpi-Si{N) for N = 0,...,L. Hence, Ai{L-N) = 
—Ai{N), and the average over a period = 1, . . . , L vanishes. □ 

Also for (3=1 and 7 = 1, the result in Example 111.1 1 agrees with Theorems 
I3.4l and 13.111 For divisor methods, we have the following general result. (By 
Example 1 1 1 . 1 1 and (|3.14p . the result does not hold for any other /3.) 

Theorem 11.3. Consider the /3-linear divisor method and suppose that /3 = 
k/2 for some integer k. Then (|3.5p holds for any rational pi, . . . ,Pm- 

Proof. By (|3.14p . the result holds for /3 if and only if it holds for /3 + 1. 
Hence it suffices to consider /3 = 1/2 and /3 = 1. The case /3 = 1/2 is part of 
Theorem HO] 

For /3 = 1, we argue similarly. (Cf. the proof of Theorem lA.llI ) Let 
< A^ < L and let t > be such that Si{N) = \pit]^, see (|2T| and (jO)) . 
Recall that [pit]^ = [pit\ except that [pit\ — 1 also is possible when pit is an 
integer. 

Since piL is an integer, if we first for simplicity assume that pit is not, 

[pii2L - t)], = 2piL - [pit]^ - 1 = 2p,L - s^iN) - 1, (11.2) 
where summing over i yields 

m m 

^[pi(2L-t)]i = 2L-^Si(A^)-m = 2L-A^-m. (11.3) 

i=l i=l 

Consequently, (1A.3P again shows that these numbers yield the seat distri- 
bution for 2L — N — m seats. (Note that L = Y17!=iPi^ — ^^'^ thus 
2L — N — m > 0.) Hence, using pi.2p again, 

Si{2L -m-N) = [pi{2L - t)]^ = 2piL - Si{N) - 1. (11.4) 

The argument works also, with a little care, if some pit is an integer, and 
shows that pi.4p holds generally, and that ties in the distribution of Si{N) 
correspond to ties in the distribution of Si{2L — m — N). Hence, we have, 
regarding Aj as a random variable when there is a tie, 

Ai(2L-m-A^) = Si{2L-m-N)-{2L-m-N)pi = mpi-Ai{N)-l. (11.5) 

Taking the expectations in case of a tie, and then the average over the 
period A^ = 0, . . . , L — 1, when 2L — m — N runs through the period L — m + 
1, . . . ,2L — m, we see that if the average is x, then x = mpi — x — 1; hence 
X = ^{mpi — 1), and thus (|3.5p holds when j3 = 1. □ 

For the 1-quota method (Droop), we have seen that (I3.16P holds in Ex- 
ample [TLTI in fact, when m = 2, it holds for any {pi,P2) by Theorem lll.3( 
because the 1-quota method equals the 1-linear divisor method when m = 2, 
see Appendix IB. II However, it does not hold for m = 3. 

Example 11.4. Let m = 3 and {pi,P2,P3) = (|)f)i)- Consider the 1- 
quota method (Droop). For A^ = 1,2,3,4,5, the smallest party gets, taking 
the average when there is a tie, 0,0, 1, i, 1 seats. (Note the non-monoticity, 
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an instance of the Alabama paradox which frequently occurs for all quota 
methods, cf. [17].) The seat excesses are — |) — §> +f > — ^7 0, with average 
— j^, while (|3.16p would give Pi ~ ^ = Thus (|3.16p does not hold. 

Note that this example shows that the heuristic argument in Remark 13.121 
is not always valid when pi, ■ ■ ■ ,Pm are rational. In this example, the re- 
tracted seat in Remark 13. 121 is not taken uniformly; a calculation shows that 
it is taken from the smallest party only ^ of the times. 

We can ask the following for any /3-linear divisor method or 7-quota 
method, where the cases 7 = 1 in Problem 111.51 and f3 = 1/2, (3=1 and 
7 = in Problem 111.61 are especially interesting. (By Theorem IA.12| for the 
linear divisor methods, it suffices to consider < /3 < 1.) 

Problem 11.5. Find a general formula for the asymptotic bias limE Aj for 
rational pi, ■ ■ ■ ,Pm- ( Theorems 1 1 1 . 2 1 and 1 1 1 . 3l yield the answers for /3 G Z/2 
and 7 = 0.) 

Problem 11.6. Find a general formula for the asymptotic distribution of 
Aj for rational pi, . . . ,Pm- 

Appendix A. Divisor methods 

A divisor method is based on a given sequence of numbers d{n), n > 1, 
with < d(l) < d{2) < d{3) < . . . . Different choices of d{n) give different 
methods. A number of divisor methods that have been used or discussed are 
shown in Table [1] The most important ones are the widely used methods by 
Jefferson/D'Hondt and Webster /Sainte-Lague. (Notation varies. In e.g. Jl|, 
the sequence is denoted d(0), d{l), . . . ; thus their d{n) is our d{n + 1).) 

Divisor methods can be described in several different ways that yield the 
same result. The methods have been invented and reinvented in different 
guises; tradition varies, and different types of formulations are used in, for 
example, election laws and other official and non-official descriptions from 
different countries. There are two main types of formulation, and we give 
each in two different versions. (The equivalence of the different types of 
formulations have been known for a long time. For example, D'Hondt pro- 
posed his method in [5] using essentially our formulation ID II below, and later 
showed the equivalence to ID4| see jisl . p. 125].) 

In the first type of formulation, the method is seen as a way of rounding. 
Given a sequence d{l),d(2), ... as above, we define the d-rounding of a real 
number x > by 

[x]^:=n if d{n) <x <d{n + l), (A.l) 

where d{0) = 0. Note that this is unambiguous only if d{n) < x < d{n + 1); 
if X = d{n) for an integer n, then both [x]^ = n and [x\^ = n — 1 are accept- 
able values. (This is important in the case of ties, see below.) In particular, 
d[n) = n — ^ (Webster's method) yields standard rounding; d{n) = n (Jeffer- 
son's method) yields rounding downwards; d(n) = n — 1 (Adams's method) 
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QcUlllldji VdlUct) I lULlllQcQ I 


Jefferson, D'Hondt 


n 


1, 2, 3, 4 ... 


Webster, Sainte-Lague 


n — h 

2 


0.5, 1.5, 2.5, 3.5, . . . 


Adams 


n — 1 


0, 1, 2, 3, ... 


Imperiali 


n+1 


2, 3, 4, 5, ... 


Danisli 


n — 1 


0.333, 1.333, 2.333, 3.333, . . . 


Adjusted Sainte-Lague 


fo.7, n = l 




1 n — |, n > 1 


0.7, 1.5, 2.5, 3.5, . . . 


Cambridge Compromise 


(n - 6)+ 


0, 0, 0, 0, 0, 0, 1, 2, 3, . . . 


Huntington 


^/n{n - 1) 


0, 1.414, 2.449, 3.464, . . . 


Dean 


2n{n-l) 
2n-l 


0, 1.333, 2.400, 3.429, . . . 


Estonia 


^0.9 


1, 1.866, 2.688, 3.482, ... 


Macau 


2n-l 


1, 2, 4, 8, ... 



Table 1. Some divisor methods. (For Dean's method, d{n) 
is the harmonic mean of n and n — 1; cf. Webster's and Hunt- 
ington's methods with the arithmetic and geometric means.) 



yields rounding upwards. (Other choices of d{n) yield non-standard round- 
ings. Note that "rounding" in general should be interpreted in a very weak 
sense, and that | [x]^ — a;| > 1 may be possible, see Remark IA.3I ) In this 
formulation, the numbers d{n) are sometimes called signposts, see |1| and 

a. 

Using the concept of d-rounding, and the notation from Section |2l the 
divisor method may be defined as follows: 

Divisor method, formulation Dl. Let 



where D is a (real) number chosen such that X^i^i ~ ^■ 

The rationale for this formulation is that D is regarded as the price of 
a seat, i.e. the number of votes that a seat "costs", which would give Vi/D 
seats to party i; this real number is d- rounded to an integer to obtain Sj. 
The price D is set by "the market" (i.e. the election officer, or computer) so 
that the desired total number of seats is distributed. This interpretation can 
be further combined with the definition (jA.ip to say that the price of n seats 
is d{n)D votes. 

Remark A.l. In particular, (lA.ip yields [x]^ = if and only if x < d{l), 
or perhaps x = d{l). Thus, if d{l) = 0, then [x]^^ > 1 for every x > 0, while 
if d{l) > 0, then [x]^ = for small positive x. 

Hence, if d{l) = 0, then (1X2]) yields Si > 1, so every party (with at least 
one vote) is guaranteed at least one seat. This is unacceptable in general 
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elections, where on the contrary there often are special threshold rules to 
prevent small parties from getting seats. However, this is acceptable, and 
may even be desirable, when distributing seats between states or constituen- 
cies, where typically there is a side condition of at least one seat each (and 
sometimes a minimum of two or more seats); a typical example is the United 
States Constitution jl|. (If d(l) = 0, we assume that N > m, so that there 
are enough seats to satisfy this minimum requirement.) 

A divisor method is by ^ called impervious if d{l) = and pervious if 
d(l) > 0. 

If we move D continuously from oo down to 0, it is seen that the total 
number of seats Yl^i given by (IA.2|) grows monotonously from (when 
d{l) > 0) or m (when d{l) = 0) to oo, and there is a value of D that yields 
the desired sum YllLi = N . In general there is an interval of such D, all 
yielding the same seat distribution Sj , and this determines all Si uniquely. In 
exceptional cases, there is only a single D that works; in this case we must 
have equalities Vi/D = d{si) for some party i and Vj/D = d{sj — 1) for some 
other party j, see (1A.2P and (lA.ip . and in this case the result (sj)™ is not 
unique; we may move one seat from party j to party i and ()A.2p still holds. 
Such cases are known as ties; they are usually resolved by drawing lots, 
but other special rules for ties are used in some countries. The formulation 
IDll thus gives a well-defined method in the sense that the numbers Si are 
uniquely determined, except possibly when there are ties. (Note that the 
number D is not uniquely determined.) 

It is sometimes, as in Section 21 convenient to rewrite ()A.2p as follows 
using X := 1/D or t := V_/D. A divisor method may thus also be called a 



multiplier method, as in jl3|. 



Divisor method, formulation D2. Let 

Si-=hx]d = [Pit]d, (A.3) 
where x or i is a (real) number chosen such that = N. 

Formulations of the type IDll of divisor methods are the standard in USA, 
where Thomas Jefferson formulated a method of this type (with rounding 
downwards) in 1792 for apportionment; Jefferson's method was adopted by 
Congress in 1792 and used until 1832 (jj. Four other divisor methods (using 
formulation IDll with different d{n), see Table [T]) have also been important 
in USA; Huntington's method is used since 1941 and Webster's method has 
been used earlier, while Adams's and Dean's methods have never been used 
but have frequently figured in discussions, see jl5| and [l|. 

In Europe, this type of formulation is unusual (although [lo| and the 



current German election law [35|, § 6 (2)] are examples), and divisor methods 
are more commonly defined by formulations such as the two following (closely 
related) ones. Note that these formulations are algorithmic, in the sense that 
they directly yield workable methods, while with IDll one has to search for 
a suitable divisor. (Preferably by computer, where such searches are easy.) 
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The European tradition of divisor methods is younger than the American 
one, see e.g. 0, 0, Hi. D'Hondt's method (which is the same as Jefferson's) 
was proposed by D'Hondt in 1882 and first adopted in Belgium 1899; it 
is now used in many countries. Sainte-Lague's method (which is the same as 
Webster's) was proposed by Sainte-Lague (H, H^l in 1910, and is also used 
in several places. 

Divisor method, formulation D3. The N seats are awarded sequentially, 
with each seat given to the party i that has the highest value of the quotient 
Vi/d{si + 1) where Sj is the number of seats that the party has received so 
far. 

The quotient Vi/d{si + 1) is called the comparative figure of the party. 
Note that the comparative figure is updated (decreased) each time the party 
receives a seat. 

Divisor method, formulation D4. Divide each Vi by the numbers d{l), 
d{2), . . . (as far as necessary). Assign the seats to the N largest quotients 
Vi/d{j). 

With formulation ID4| the quotients Vi/d{j) can be arranged in a table 
(matrix); alternatively, the sequences of quotients obtained for different par- 
ties can be merged into a single sequence in decreasing order. It is immediate 
that formulations ID3 l and ID4] are equivalent; ID3I picks the largest quotients 
in lD4l in decreasing order. 

Of course, ties can occur in formulations ID3l and lD4l as well, and again they 
are resolved by drawing lots, or possibly by some other rule. (In formulation 
ID3| ties may also occur at intermediate stages, but they do not affect the 
final result.) 

To see that the formulations ID1HD4I are equivalent and yield the same 
method, note first that (1A.2P can be written, using the definition (lA.ip . 

d{si)<^<d{si + l), i = l,...,N, (A.4) 

or, equivalently (interpreting Vi/0 = +oo) 

^u^,<D<-^ i = l,...,N. (A.5) 
d[Si + l) d{Si) 

Consequently, (sj)^]^ are such that 

max— — < mm— — -; (A. 6) 

i d{si + l) - t d{si)' ^ ' 

conversely, for any such Si we may choose D such that (lA.Sp holds. Thus 
formulation IDll is equivalent to: 

Divisor method, formulation D5. si, . . . , Sm are chosen such that 
YX=i Si = N and dAH) holds. 



Furthermore, it is easily seen that formulation ID5l is equivalent to lD4| and 
thus also to ID3I Hence all formulations are equivalent. (It is easily verified 
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that also ties appear simultaneously in the different formulations, except 
that in formulation ID3I there may be additional intermediate ties that do 
not affect the result.) Further equivalent formulations are given by [s]. 

We see also, from (IA.5|) . that the possible choices of D in formulation ID II 
that yield si, . . . ,Sm are D S where 

D_ := max , D+ := min (A. 7) 

I d[Si + 1) t d{Si) 

Thus, -D+ equals the winning comparative figure when the last seat is awarded 
in formulation ID3| and equals the winning comparative figure for the next 
seat (if there were one) . 

Remark A. 2. It follows immediately from any of the formulations ID1^D5I 
that the sequences d(n) and cd{n), for any constant c > 0, define the same 
divisor method. For example, Sainte-Lague's method is given in Table [1] 
with d{n) = n — ^; it can as well be defined with d{n) = 2n — 1. (In fact, 
the method is traditionally defined with this sequence, see e.g. Sainte-Lague 
2^, l^l- It is therefore also called the odd-number method.) Similarly, the ad- 



justed Sainte-Lague's method (used in Sweden and Norway) is traditionally 
defined with the sequence 1.4, 3, 5, 7, ... , and the Danish method with the 
sequence 1, 4, 7, 10, ... . (This method is used in Danish parliamentary elec- 
tions for part of the distribution of adjustment seats between consituencies; 
the main distribution among parties is by the method of largest remainder.) 
It may be convenient to use integers; on the other hand, Imperiali's method 



(used in Belgian local elections) is in the election law |33l . Art. 56] described 
using 1, li, 2, 2i, . . . , i.e., {n + l)/2. In Table H we have normalised d{n) 
to the form n -|- /? whenever possible, cf. Section lA.ll 

Remark A. 3. Several papers impose the condition that n — 1 < d{n) < n, 
but this is really not necessary. This condition is natural for the interpre- 
tation of the method as a kind of rounding in formulation IDH since it is 
equivalent to [n]^ = n for any integer n S N, see (jA.ip . In particular, this 
means that if the exact proportions qi = piN happen to be integers, then the 
divisor method gives Sj = Qi. (Take D = V/N in (IA.2P ; cf. (|2.ip .) However, 
the method is well-defined also for other sequences d{n), and as seen in Ta- 
ble [H there are some divisor methods currently in use for general elections 
that do not satisfy this condition (the Imperiali method [sll. Art. 56] and the 



methods in Estonia |34l . § 62 (5)] and Macau [36, Artigo 17]). Formulation 
ID II remains formally valid in these cases too, but in practice these methods 
are described using e.g. formulation ID3I where there are no problems of in- 
terpretation for any increasing sequence d{n). (Methods satisfying Si = qi 
when all q^ are integers are called weakly proportional by [l|.) 

Remark A. 4. We required above that the sequence d{n) is strictly increas- 
ing. In fact, it can be defined assuming only < d{l) < d{2) < . . . , with 
only a few minor modifications above. This extension is hardly used in prac- 
tice, except that a minimum requirement of at least r seats for each party 
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for some r > 1 can be achieved by taking d{l) = d{2) = • • • = d{r) = 0, cf. 
Remark lA.ll for the case r = 1. (We then need A'^ > rm.) 

For example, the Cambridge Compromise (lo| proposed for apportion- 
ment of the European Parhament can be described as giving each state a 
"base" of 5 seats and distributing the rest by Adams's method (subject to a 
maximum of 96 seats for each state); this is (e.g. by formulation ID1| see also 
Theorem I A . 1 2 1 below) equivalent to giving each state 6 seats and distributing 
the rest by D'Hondt's method (subject to the same maximum). If we ignore 
the restriction to at most 96 seats (as done in Table [T] for simplicity), this is 
the same as a divisor method with d{n) = for n < 6 and d{n) = n — 6 for 
n > 6, i.e. d{n) = [n — 6)+. The maximum restriction is easily implemented 
if e.g. formulation ID3I is used, cf. (lo| : formally we can also implement it by 
redefining d{n) = oo for n > 96. 

Remark A. 5. The name "divisor method" is used for methods of this type 
in any of the formulations above. However, the word divisor is usually used 
for the (variable) number D in formulation ID1[ but for the (fixed) numbers 
d(l), d(2), ... in formulations [D3tiD4l 

A.l. Linear divisor methods. We say that a divisor method is linear if 
d{n) = an + h for some a > and 6 G M. By Remark IA.2[ we may replace 
d{n) by d{n)/a = n + b/a. We define (3 := b/a + 1. In other words, we may 
assume that a = 1 and 

d{n)=n-l + /3 (A.8) 

for some real /3. The reason for our choice of f3 as parameter, and thus 
formula (lO]) . is ^^^ below. Note that f3 = d{l). 

We call the divisor method with d{n) given by (IA.8|) the /3-linear divisor 
method. The par ameter /3 is called proportionality index [Proportionali- 
tdtsindex) by 1^. The method is called q-stationary multiplier method 
(where q = f3) by [13 1. 



With a linear divisor method (lA.Sp . the definition (jA.ip of d- rounding 
yields that if [x]^ = n, then n — l + /3<x<n + 13, oi equivalently x — (3 < 
n < X — (3 + 1. Consequently, d-rounding equals the /3-rounding defined in 
and (|23]); i.e. 

= N/3 , (A.9) 

at least provided x > and x > /3 — 1. Hence, the /3-linear divisor method 
can, using formulation IDH be defined by (12. 7p as we did in Section [2] If 
f3 > 1, we here have to assume that A^ is so large that Vi/D > (3 — 1 and 
thus [vi/D]g > 0. 



Remark A. 6. If /3 > 1 and A^ is small, (12. 7p may give a negative seat 
number Si for a very small party. Of course, this has to be replaced by 0, 
as given by (lA.2p . Hence, we may still use (12. 7p if we ignore any party that 
would get a negative number of seats. 
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If [vi/D\Q < —1 for some i, then Vi/ D < P — 1 hy (12. 5p and thus for any- 



party j, using (1231), 



D 



< 



/9 



-iP-l) 



< 



Pi 

Pi 



Hence, if > 



and summing over j we get N < {1/pi — m){f3 — 1 
(1/minjpj — m){f3 — 1), then (|2.7p holds without exception. 

Remark A. 7. The definition above requires /3 = d{l) > 0. We can extend 
it to /3 < by replacing (jA.Sp by 



d{n) = (n - 1 + /3)_ 



(A.IO) 



i.e., letting d(n) = for n < 1 + [\(3\\, cf. Remark |X1 Then (1X9]) still 
holds for X > 0, and (12. 7p holds as before. By Remark IA.4[ every party now 
gets at least [|/3|J + 1 seats, and thus we need A^ > m{[\(3\ + Ij). 

The Cambridge Compromise (without maximum rule) is an example, with 
(3 = -5. 

Remark A. 8. Linear divisor methods are often considered only for the case 
< /3 < 1, which is equivalent to the condition n — 1 < d{n) < n discussed 
in Remark IA.3I However, we see that they are well-defined for arbitrary real 
(3, although some care may be needed for small A^. Note that at least the 
case /3 = 2 is used in practice. 

In the main part of the paper we consider only linear divisor methods. 
Note that many of the methods in Table [1] are linear: Jefferson/ D'Hondt 
(/S = 1); Webster/Sainte-Lague (/3 = 1/2); Adams (/3 = 0); Imperiah {/3 = 
2); Danish (/3 = 1/3); Cambridge Compromise (/3 = —5). 

Remark A. 9. The adjusted Sainte-Lague method also has a linear d{n) as 
in (lA.SP (/3 = 1/2) except for d(l), which does not affect asymptotic results. 
More precisely, the adjusted Sainte-Lague method differs from Sainte-Lague's 
method only in the value of d(l), and [x]^ is the same for both methods for 
all X > 0.7; hence they give exactly the same result as soon as every party 
gets at least one seat by the adjusted method. (But the adjustment makes it 
more difficult for a small party to get the ffist seat.) For asymptotic results 
as in the present paper, there is thus no difference between the adjusted 
Sainte-Lague method and Sainte-Lague's method. 

Remark A.IO. Huntington's and Dean's methods are asymptotically linear 
in the sense that d{n) = n— l+j3+o{l) as n — t- oo, in both cases with /3 = 1 /2. 
An argument similar to the proof in Section |4] shows that as A^ — > oo, the 
probability that these methods yield the same result as Webster's tends to 
1. In particular. Theorem 13.71 (with /3 = 1/2) holds for these methods too, 
and Theorem 13.41 shows that they are asymptotically unbiased. 
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Note that our asymptotic approach does not answer the controversial ques- 
tion whether Webster's or Huntington's method is the most fair and unbi- 
ased, see [H; this depends mainly on how one measures the bias for small 
parties (states), in particular for the ones obtaining just 1 seat. 

Although linear divisor methods may be biased, see Section |3l they are 
perfectly proportional with respect to changes in the total number of seats. 
(If /3 < or /3 > 1, we assume that N is not too small, see Remarks EIHAS) 

Theorem A.ll. Consider a 13 -linear divisor method. If the number of seats 
is increased from N to N + L, and Lpi, . . . , Lpm all are integers, then party 
i gets exactly Lpi seats more: 

Si{N + L) = Si{N) + Lpi. (A.ll) 

Proof. This is easiest seen using the multiplier formulation ID2I If t yields 
Si{N) = [pit]^ with YT=i ^ii^) = tlien ^ + ^ yiel'is 

Si = [Pi{t + L)]p = \pit\i^ +PiL = Si{N) +PiL, 
with J:T=i s^ = YZi s^iN) + TZlP^L = N + L. □ 

There is also a simple relation between the methods for two different values 
of (3 that differ by an integer. 

Theorem A. 12. Consider the (3-linear divisor method and suppose that N 
is so large that every party gets at least one seat. ( Otherwise, ignore the 
remaining parties.) Then the (3-linear divisor method yields the same result 
as first giving one seat to each party and then distributiong the remaining 
N — m seats by the [(3 + l)-linear divisor method. 

Proof. This follows immediately from formulation ID4| noting that by (jA.SP , 
or more generally (jA.lOp . dpj^i{n) = dp{n + 1) (using /3 as a subscript to 
denote the divisor sequences). 

Alternatively, this follows from (|2.7p and the relation [a;]^_|_]^ = [x]^ — 1, 
which follows from (|2.6p . □ 

By repeated applications of this theorem, we can reduce any /3-linear di- 
visor method to the case < /3 < 1. 

Example A. 13. Adams's method {(3 = 0) is the same as giving each party 
one seat and distributing the rest by Jefferson's method {(3 = 1). 

Example A. 14. The Imperiali method {(3 = 2) is equivalent to using 
D'Hondt's method (/3 = 1) but retracting one seat from each party and 
redistribute them (among the parties that had at least one seat) by contin- 
ued application of D'Hondt's method. 

Example A. 15. The Cambridge Compromise (/3 = —5), is as said in Re- 
mark IA.4I the same as giving each party (state) 5 seats and distributing the 
rest by Adams's method (/3 = 0), or as giving each party (state) 6 seats and 
distributing the rest by Jefferson/D'Hondt's method {(3 = 1). 
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A. 2. Divisor methods and proportionality. We make a formal defini- 
tion of proportionality for election methods, using an asymptotic property, 
since no election method is exactly proportional except in exceptional cases. 
(Balinski and Young uses a different definition, cf. Remark IA.3I ) 

Definition A. 16. An election method is (asymptotically) proportional if, 
for any number m of parties and any proportions of votes pi, ■ ■ ■ ,Pm, for 
each party i, 

Si{N)/N^Pi as N ^ oo. (A. 12) 

Many proportional methods that are used are divisor methods, but not all 
divisor methods are proportional. In order to characterise the proportional 
divisor methods, recall that a positive measurable function / on some interval 
{A, oo) is regularly varying with index p, where p is a real number, if 

f{\x)/f{x)^\P asx^oo, (A.13) 

for every A > 0. In the special case p = 0, i.e. when f{\x)/f{x) — t- 1 as 
j; — )• oo for every A > 0, we say that / is slowly varying. A typical example 
of a slowly varying function is logj;. A sequence c„ is said to be regularly 
varying with index p (and slowly varying if p = 0) if the function cui is; this 
is equivalent to c\xn\/(^n — )• A'' as n — )• oo, for every A > 0. (See [3[, where 
many more results are given.) 

Theorem A. 17. For a divisor method defined by d{l),d{2), . . . , the follow- 
ing are equivalent. 

(i) The method is proportional. 

(ii) For any m and pi, . . . ,Pm, o-nd all i,j < m, 

^A^^El asN^oo. (A.14) 

s^iN) p, ^ 

(iii) x I— )■ [x]^ is regularly varying with index 1. 

(iv) X I—)- [xj^/x is slowly varying. 

(v) The sequence d{n) is regularly varying with index 1. 

(vi) The sequence d{n)/n is slowly varying. 



The proof shows that it suffices that (IA.12p . or (ii) holds when m = 2. 



Proof. [(I)] <S=^ [(ii)] If (IA.12|) holds for each i, then also (IA.14|) holds. On 
the other hand, if (IA.14[) holds, then summing (IA.14|) over all j we obtain 
N/si{N) -^l/pi, and thus flAH]) . 



m 



By (IA.2[) . Si = [vi/D]^ = \pix\^ with x := V/D, and hence 
can be written 

h^^Pl asx^oo, (A.15) 

since = Yll^i — YliLi \Pi^]d ^^"^ tl^us N ^ oo <^=^ x — )• oo. If 
[x]^ is regularly varying with index 1, then (lA.lSP follows from (IA.13P with 
A = Pj/pi- Conversely, if (ii) holds, and A > 0, consider the case m = 2 and 
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pi = \/{l + X),P2 = 1/(1 + A). Then ( jAlSl) with = (1,2) and y:= p2X 
yields [Xy]^ / [?/]cZ ~^ A, which is (IA.13|) with p = 1- 
(iv)[ Directly by the definition (IA.13|) . 

Let g{x) be the inverse of [x]^ defined by g(x) := inf{y : 



111 



111 



[y]^ > x}. Then (IA.1[) shows that g{x) = d{[x\ + 1). jsl, Theorem 1.5.12] 
shows that if [x]^ is regularly varying with index 1, then so is g{x). This 
implies d{n + l)/d{n) — )• 1 as n — )• oo, and it follows that d{n) is regularly 
varying with index 1. The converse follows in the same way, now using the 
inverse of d([a;J) which is [x]^ + 1. 

Directly by the definition. □ 



(v) 







VI 



In particular. Theorem lA. 171 shows that any divisor method with 

d(n)/n— )-a as n — )• oo, (A. 16) 

for some a > 0, is proportional. Of the methods in Table [H the only methods 
that do not satisfy (IA.16|) . and thus are proportional, are the methods of 
Estonia and Macau. In fact, it is easily seen from formulation ID3l or ID4l that 



the Estonian method is the same as D'Hondt's method applied to v^^'^ = 
^1.111.... i^jjyg makes Sj proportional to -y^m---. This gives a small but clear 
advantage to larger parties. Macau's method, on the contrary, favours small 
parties and encourages splintering. (With only 12 directly elected members 
of the Legislative Assembly of Macau, it essentially imposes a maximum of 
2 seats per party.) 

In fact, the proportional divisor methods in Table [1] and all linear divisor 
methods, satisfy the following stronger form of proportionality. 

Theorem A. 18. For a divisor method with d{n) = an + 0(1) for some 
a > 0, Si{N) = piN + 0(1). (The implicit constant may depend on m, hut 
not on anything else.) 

Proof. By Remark IA.2[ we may replace d[n) by d[n) /a, and thus assume 
d{n) = n + 0{l). 

If d{n) = n + 0(1), then also d{n — 1) = n + 0(1) and thus (lA.ip shows 
that if [x]^ = n, then x = n + 0(1); in other words, 

[x], = x + 0(l). (A.17) 

Consequently, (1A.2P yields, with y := V/D, 

Si = ^ + 0{l)=piy + 0{l). (A.18) 

Summing over i we find = 'YllLiPiV + ^(1) = ^ + ^(1), and thus (jA.lSP 
yields Si =piA^ + 0(1). □ 

A. 3. Uniformity. Divisor methods satisfy the following consistency prop- 
erty, called uniformity in [ij, when we consider subsets of the set of all 
parties. This is obvious from any of the formulations ID1HD5[ and would 
perhaps not be worth mentioning, except for the fact that it does not hold 
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for quota methods; in fact, under some weak extra conditions, it is satisfied 
only by divisor methods jlj. 

Theorem A. 19. Suppose that some set of parties Ai, . . . , Ai (where 1 < i < 
m ) get together N' seats in an election by a divisor method. Then the number 
of seats for each of these parties is the same as if N' seats were distributed, by 
the same divisor method, in an election with only these parties participating 
(and obtaining the same numbers of votes Vi). □ 



Appendix B. Quota methods 

In a quota method one first determines a quota Q; different quota methods 
differ in the formula for Q, see below. {Q can be seen as the standard price 
of a seat, just as the divisor D in the formulation IDll of divisor methods.) 
The seats are then distributed as follows: 

Quota method, formulation Ql. Divide the numbers of votes Vi by Q, 
and give first each party as many seats as the integer part [vi/Q\ of its 
fraction. Any remaining seats are given to the parties that have largest 
fractional part {vi/Q} = Vi/Q— [vi/Q\. (Or, equivalently, largest remainder 
at the division Vi/Q, since the remainder is Vi — [vi/Q\Q, which is Q times 
the fractional part Vi/Q — [vi/Q\.) 

The most common quota method is the method of largest remainder, or 
Hamilton's method, in Europe often called Hare's method, which uses the 
simple quota, also called Hare quota, Q := V/N , i.e. the average number 
of votes per seat. Note that then Vi/Q = Qi, the exactly proportional real 
allocation, see (12. ip . (Alexander Hamilton suggested the method in 1792 
for apportionment to the US Congress; it was approved by Congress but 
vetoed by president Washington; it made a comeback and was used (under 
the name Vinton's method) 1850-1900, see [l|. The namn Hare's method 
is also well-established but is a misnomer; Thomas Hare really advocated a 
different method, the Single Transferable Vote, STV, which is not based on 
party lists and is not treated here.) 

Another version is Droop's method which uses the Droop quota Q := 
V/{N + 1). Also the Imperiali quota Q := V/{N + 2) has occasionally been 
used (e.g. in Italy until 1993 ^]). 

We define, more generally, for any real number 7, the ■j-quota method to 
be the quota method with quota 

(For 7 < 0, we assume that N > \'y\.) Thus Hamilton/Hare's method 
is the case 7 = 0, Droop's method is 7 = 1, and the Imperiali quota is 
7 = 2. (Kopfermann fl8| calls the 7-quota method Rundungsverfahren mit 
dem Proportionalitdtsindex p, where his /f = (7 + l)/2, i.e., 7 = 2/9 — 1.) 
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Remark B.l. The quotas above are sometimes rounded to integers (up- 



wards, downwards or by standard rounding); see |22| for several different 



examples in current use. There is no mathematical reason to round the 



quota; on the contrary, there are reasons against. Kopfermann |18| points 
out that quota methods involve rounding at the end, but that, as a gen- 
eral principle, intermediate values in a calculation should not be rounded. 
A specific problem is that rounding the quota means that the method is 
no longer homogeneous (see Section [2]). Moreover, with a rounded quota, 
in some cases a party may lose a seat by getting an additional vote, which 
ought to be unacceptable; see jl6| for a simple example with Droop's method 
(although there formulated for STV, which in the example gives the same 
result). In a general election with a large number of votes, and therefore a 
large quota, rounding the quota will usually make no difference, but in the 
exceptional cases when it does make a difference, it can be harmful and lead 
to undesirable results. 

In the present paper, we consider the methods with unrounded quotas. 



Remark B.2. In formulation Ql it is implicitly assumed that N — m < 
Yl'i^il^i/Ql — t^^t number of remaining seats is non-negative 

and not larger than m, the number of parties. For the 7-quota method, this 
is always satisfied if — 1 < 7 < 1 (including Hamilton/Hare), and also if 
7 = 1 (Droop) except in a very special case when all parties tie for the last 



seat. For 7 < —1 or 7 > 1, formulation Ql may have to be amended; in 
such cases we interpret it by using one of the equivalent formulations Q2 - Q4 
below. (These formulations are well-defined for all Q and N, up to the usual 
possibility of ties.) 



The rule in |Ql| that the remaining seats are distributed to the parties with 
largest remainders Vi/Q—[vi/ Q\ can also be expressed by saying that we find 
a suitable a and round Vi/Q upwards if the fractional part is greater than 
a, and downwards if the fractional part is smaller. In other words, recalling 
the notion of a-rounding in Section 121 a quota method can equivalently be 
described as follows. 



Quota method, formulation Q2. 

^""^ , (B.2) 



where a is chosen such that X^i^i ~ ^■ 

For the 7-quota method we have Q = V/{N + 7), and thus Vi/Q = 
ViV IQ = Vii]^ + 7)- Hence we have also the following formula. 

Quota method, formulation Q3. The 7-quota method is given by 

Si ■= \{N + 7)p4 , (B.3) 

where a is chosen such that ^Y^=\ — 
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Formulations Q2 - Q3 are the same as (12. Sp . 

Quota methods can also be described using a sequential allocation of seats 
by comparative figures as in formulation ID3l of divisor methods. (This is not 



the usual way to define quota methods, but it appears in proofs in |27[ and 



23j.) It is easy to see that Q2 is equivalent to the following. 



Quota method, formulation Q4. The N seats are awarded sequentially, 
with each seat given to the party i that has the highest value of the difference 
Vi — SiQ where Sj is the number of seats that the party has received so far. 

It is easy to see that every 7-quota method satisfies Si{N) = Npi + 0(1), 
uniformly in all m and pi, . . . ,pm, and hence is proportional, cf. Defini- 
tion IA.16I and Theorem IA.18I We omit the easy proof since we show a 
precise version of this in Theorem 13.101 

Moreover, a 7-quota method is perfectly proportional with respect to 
changes in the total number of seats, just as a /3-linear divisor method, 
cf. Theorem lA.llI 

Theorem B.3. Consider a 'j-quota method. If the number of seats is in- 
creased from N to N + L, and Lpi, . . . , Lpm all are integers, then party i 
gets exactly Lpi seats more: 

Si{N + L) = Si{N) + Lpi. (B.4) 

Proof. Obvious from ()B.3p . since the same a works for N and N + L. □ 

B.l. Two parties. In the simple case m = 2, when there are only two 
parties, the /3-linear divisor method and the 7-quota method with 7 = 2/3 — 



1 coincide |l8l . Satz 6.2.6]. (In particular, then Webster /Sainte-Lague = 
Hamilton/Hare and Jefferson/ D'Hondt ~ Droop.) In fact, it is easy to see 
that 

:= [(iV + 2/3 - , i = l,2, (B.5) 

gives si -|- S2 = A^, (at least for some choice in case of a tie). Thus, for the 
/3-linear divisor method we can take D' = 1/{N + 2/3 — 1) = 1/{N + 7) and 
D = D'V = V/{N 2/3 - 1) in (iTT]) . and for the 7-quota method we can 
take a = 13 in (12. Sp . Hence, both methods yield the result given by (IB.Sp . 

Appendix C. Lemmas on uniform distribution 

Proof of Lemma \5.1\ First, by Lemma k.lj the sequence 

{{nyi/m + ai/m}, {nyk/m + ak/m}) (C.l) 

is uniformly distributed in [0, 1)'^. Thus, by the definition of Modm, 

{Modm{nyi + ai), . . . , Modm{nyk + a^)) (C.2) 

is uniformly distributed in [0, m)'^. If a sequence (x„) is uniformly distributed 
in [0,m), then the sequence of pairs ({x„}, [xn\) is uniformly distributed in 
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[0,l)x{0,...,m — 1}. This extends to k dimensions, and thus it follows 
from (IC.2P that if £nj '■= Modm(L^yj + "^jJ) = [^odm{nyj + aj)\, then 

{{nyi + ai}, . . . ,{nyk + ak}Jni, ■ ■ ■ Ank) (C.3) 

is uniformly distributed in [O,!)*^ x{0,...,m— 1}'^. 

The argument just given shows that (IC.3P is uniformly distributed also 
along any subsequence n = mu + no, f > 1. We have = Vio6.m{n — 
X^j=i^nj)) and it follows that along any such subsequence, ^„ is determined 

by ini, ■ ■ ■ ,£nk- Consequently, the vectors (IS.lip are uniformly distributed 
along any such subsequence, and therefore also along the full sequence. □ 

The following lemma is essentially taken from (l2| . (There Riemann in- 
tegrability is assumed, but as the proof below shows, Lebesgue integrability 
suffices.) 

Lemma C.l. Let X = {Xi, . . . ,Xk) be a random variable with an absolutely 
continuous distribution in R'^. Let Vn be any sequence of constants with 
— 7- oo. Then 

{{VnXl],...,{VnXk],Xr,...,Xk) A (C/i,...,C/fc,Xi,...,Xfc) (C.4) 

where Ui, . . . ,Uk ~ U(0, 1) are independent of each other and of X. 

Proof. The assumption that X is absolutely continuous means that it has a 
density function /, which is a Lebesgue integrable function on M'^. 

We regard {unXj} and Uj as elements of the circle group T = M/Z. Thus 
the random vectors in (1C.4P are elements of the group T'= X M'^, and by 
standard Fourier analysis, it suffices to show that for any integers ii, . . . ,ik 
and real numbers ti, . . . ,tk, 

Since i^nXj — {i^nXj} is an integer, we have, with / as above, 

= EgJ^. i(2-f..n+t,)x, ^ /(27r4z.„ + ti, . . . , 27r4z^n + tfc). (C.6) 

If (£i,...,4) = (0, ...,0), then (|U31) is trivial. If (4,..., 4) / (0, ...,0), 
then |(27r£iz^„ + ti, . . . , 27r£fcZ^„ + tfc)| — )• oo, and thus (IC.6P and the Riemann- 
Lebesgue lemma show that the left-hand side of (IC.SP tends to 0, which 
verifies (IC.SP in this case too, since 

because E e^'^^^^^^ = when / 0. □ 
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Appendix D. Moments of order statistics for random pi 

Let Pi, . . . ,Prra be random and uniformly distributed on as in Section |6l 
and let < • • • < P(m) be the order statistics. We give here, for complete- 
ness, a calculation of moments of by a well-known method. Note that we 
here follow standard convention for order statistics in probability theory and 
order pi in increasing order, in contrast to Section [6] where we order them in 
decreasing order; thus pj^j = P[m+i-k)i which should be remembered when 
using the results below in Section [6l 

Let Ti , . . . , Tra be m independent identically distributed random variables 
with the exponential distribution Exp(l), and let Z := X^I^i Then 
{Ti/ Z, . . . ,Tm/ Z) is uniformly distributed on the simplex ©m, and more- 
over independent of Z. This means that we can take pi = Ti/Z. For the 
corresponding order statistics we thus have p(^) = Tq,-^/Z] moreover p(^) and 
Z are independent and thus, for any k and 



so the moments of pQ^-^ are given by 



Furthermore, Z := Xll^i ^« bas the Gamma distribution T{m) with moments 



EZ^= ^^Z^f^ =m(m + l)---(m + £-l). (D.3) 
Tim) 



To find the moments of the order statistics T^^-j, we regard Ti, . . . , Tm as 
the times of the first events in m independent Poisson processes. (Alterna- 
tively, we may regard them as life-lengths of m identical radioactive atoms.) 
It is then well-known, and easy to see by the lack of memory in exponen- 
tial distributions, that the increments (or waiting times) Vj := T(^j^ — Tq-i) 
(with T(-Q) := 0) are independent exponential random variables with Vj ~ 

Exp(l/(m -|- 1 — j))- Hence T(fc) = Yl^=i with these Vj, and moments are 
easily computed. In particular. 



1 



m + 1 — j 



m -. 

E 7- 



(D.4) 



i=-m— fc+1 
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(This can be written as — Hm-k, where Hm = X^I^i ^/^ is the m:th 
harmonic number.) Similarly, 



k Ik 

v2 



k ^ ( ^ \ 

~ ^ (m + 1 - j)2 ( ^ m + 1 



i=m— fc+l \i=m— fc+1 



Hence, by (ICT-flCT. 



i=m—k+l 

2_ ^^(1) _ 1 1 ^ 1 / ^ iV 

^ ' m(m + 1) m(m + 1) . ^-^ i-^ m{m + 1) I ^—^ i j ' 

(D.7) 

Var(p(fe))=E(p(fc))2-(Ep(fc))2 

1 ^ 1 1 / ^ iV 

m(m + 1) f2 m^fm+l) I ^-^ i j 

^ ' i=m-k+l ^ ' \i=m-k+l / 

Covariances are computed similarly. U 1 < k < i < m, then 

k k m ^ 

Cov(r(,),r(,)) = EVarl-,=E7;— Y— W= E -2 

j=l j=l ^ •' ' i=rn-k+l 

and, calculating as in (ID.SP and (ID.Gp . 



Cov(p(fc),p(^)) 

= / , y ^-^-^ — r V - V -. (D.io) 

m(m + l) ''-^ m^(m + l) i i 

^ ' i=m-k+l ^ ' i=m-k+l i=m-l+l 

Example D.l. For m = 3, the covariance matrix of (p[i],P[2])P[3]) = 

(P(3),P(2),P(1)) is 



1 
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13 -8 -5^ 
-8 7 1 
-5 1 4, 
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